FERMIONIC STRUCTURE IN THE SINE-GORDON MODEL: 
FORM FACTORS AND NULL- VECTORS 
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Abstract. The form factor bootstrap in integrable quantum field theory allows 
'' one to capture local fields in terms of infinite sequences of Laurent polynomials 

called 'towers'. For the sine-Gordon model, towers are systematically described by 
fermions introduced some time ago by Babclon, Bernard and Smirnov. Recently 
the authors developed a new method for evaluating one-point functions of descen- 
dant fields, using yet another fermions which act on the space of local fields. The 
goal of this paper is to establish that these two fermions are one and the same 
' object. This opens up a way for answering the longstanding question about how 

to identify precisely towers and local fields. 

1. Introduction 
The famous sine-Gordon (sG) model is described by the action 
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In this paper we use the parameter 



following the convention in our previous works m [21 [3] . 

In our opinion, the sG model is an ideal playground for developing new meth- 
ods of Integrable Quantum Field Theory (IQFT). On the one hand, this model is 
sufficiently complicated. Its spectrum contains, together with usual particles and 
their bound states (breathers), topologically non-trivial particles: solitons. On the 
other hand, the model is intimately related with the simplest non-trivial quantum 
affine algebra Uq{sl2), so, its study does not involve purely technical complications 
coming from considering quantum groups of higher rank. The latter property is 
closely related to the fact that in the classical case the (quasi)-periodic solutions 
of the sG equation are related to hyper-elliptic Riemann surfaces, which represent 
the simplest, but still non-trivial example of algebraic curves. All that should be 
clearly understood by a researcher who is interested in serious investigation of the 
sG model. 

During the period of rapid development of IQFT, several important results were 
obtained for the sG model. They include the discovery of the exact S-matrix [1], 
formulation of axioms for the form factors and solving equations constituting these 
axioms El El- Around the same period there appeared the work [5] which is 
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very important for us. This paper investigated the equivalence of the form factor 
bootstrap to the Operator Product Expansion (OPE) appearing in the context of 
Perturbed Conformal Field Theory (PCFT) combined with the knowledge of the 
one-point functions. 

Let us discuss the form factors in the sG model. In the original works [5l [7] they 
were constructed for operators most relevant to physics: energy-momentum tensor, 
topological current and disorder operators. With the rapid development of Confor- 
mal Field Theory (GET) which started with the famous paper [TU], the following 
interesting question arose: Find the form factors of all the local operators which 
are described in the ultraviolet limit by GET. The papers [HI [T2] were important 
in understanding this problem together with the paper |H] cited before. Due to the 
latter paper it became clear that the ideal object for this study is not the space 
of descendants of the degenerate fields, but rather the space of descendants of the 
generic primary field 



The point is that considering generic v and a we avoid resonances, and the correspon- 
dence between the sG operators and their ultra-violet GET counterparts becomes 
one-to-one. According to this logic, the first task is to compute the form factors 
of the exponential operators $a(0). This was done in the paper [9]. This paper 
required the knowledge of zero-particle form factors (one-point functions) in infinite 
volume found in [13]. It should be said also that the method of [9] originates in the 
algebraic study of correlation functions and form factors for lattice models [H]. 

It has been said that the sG model possesses the advantage of being sufficiently 
complicated while avoiding unnecessary difficulties of purely algebraic nature. Here 
we want to explain what the words "sufficiently complicated" mean. There are mod- 
els of IQFT with much simpler, diagonal, S-matrices (sinh-Gordon model, Lee- Yang 
model, model, and many others). For these models the form factor bootstrap 
considerably simplifies. However, to our mind it simplifies too much making obscure 
the mathematical structure of the solution. When the answer is simple it may allow 
different accidental representations. The sG model is sufficiently rigid for form fac- 
tors. Irrespectively of the method used to derive them, one obtains essentially the 
same formulae given by a certain integral transformation. The integrals involved in 
this transformation can be understood as quantum deformation of the hyper-elliptic 
Abelian integrals [15] . The analogy with the classical case was crucial for realising 
this fact. While quite useful, this analogy is not so straightforward. In classical 
mathematics, when the solution to a differential equation is given by an integral 
transformation, different solutions are parametrised by different contours. After 
quantisation the differential equations are replaced by a certain Riemann-Hilbert 
problem (see Section [2]), and the contours are replaced by polynomials which one 
can insert under the integral. These polynomials label different local operators. In 
this paper we consider only soliton form factors for the local operators which do not 
change the topological charge. So, the number of particles is even, say In. More 
precisely, with each n there is associated a polynomial entering the integral formula 




is introduced for convenience. 



FERMIONIC STRUCTURE IN SINE-GORDON MODEL 



3 



for 2n-particle form factors, and these polynomials are mutually related by a certain 
recurrence relation. Hence a local operator is represented by an infinite sequence of 
such polynomials which we call a tower. The precise definition is given in Section HI 

We would like to make one more historical comment. When the method of count- 
ing local fields became clear, the following important problem arose. It is well known 
that for special values of a, which correspond to degenerate primary fields, some 
descendants vanish. Conventionally these vanishing descendants are called null vec- 
tors. On the other hand the number of towers corresponding to descendants of a 
primary field is independent of a. The only solution to this apparent contradiction 
may consist in vanishing of the integrals defining the form factors in some cases. 
This is indeed the case. There are different reasons for integrals to vanish, the most 
important being the Riemann bilinear identity for quantum Abelian integrals [T6] . 

Returning to the general case, let us fix a generic a and explain how the space of 
towers is organised. A convenient language for that was introduced in the paper [T7j . 
We slightly generalise the results of this paper and change the notation. The local 
integrals of motion act on the local operators by commutators. The form factors 
obtained by this action are easy to compute, so, we shall ignore these descendants 
by local integrals of motion. The real problem is to describe the quotient space. To 
this end, following [T7] we introduce the fermionic creation operators V^2j-i' X2j-i' 
Xlj-i acting on towers, j = 1, 2, ■ ■ • . The operators 'ip2j-iJ ^2j-i correspond 
to the right chirality, and the operators ip2j-ii X2j-i ^'^ ^^e left chirality. We call 
these fermions Babelon-Bernard-Smirnov (BBS) fermions. 

Let Mq*^ = {M'^'^^jJ^Q denote the tower corresponding to the primary field 
The space of towers L^^ corresponding to local operators Oq, is obtained acting on 
Mq*'* by integrals of motion and BBS fermions. The latter must satisfy the restriction 

(1.3) #(^*) + #(V^*) = #(X*) + #(X*) • 

Later we shall say that the towers satisfying flL3p have charge 0. It is also useful to 
introduce the weight of a local operator Oq, by 

m = UMr) - Mx*) + #(x*) - Mr)) ■ 

The local operators of weight m correspond to ^^,2m,— ^ Virasoro descen- 

dants. 

As we mentioned already, there is a one-to-one correspondence between the op- 
erators in the sG model and in the corresponding ultraviolet CFT. There is certain 
arbitrariness in choosing the latter. We prefer to split 2 cos(/3£^(x)) in f ll.ip into the 
sum e~^^'^^^^ + e*^'^'-^'^ giving the first term to the CFT action and considering the 
second one as the perturbation. The CFT in question is nothing but the complex 
Liouville model. The fields ^^izji are degenerate. For m > they correspond to the 
fields from the first row of the Kac table. So, certain descendants of these fields must 
vanish. This circumstance was the main subject of [17]. Namely, it was shown that 
the Riemann bilinear identity and some additional simpler properties of quantum 
Abelian integrals imply certain relations between the fermionic descendants for the 
degenerate fields. This will be discussed in Section [HI Taking these null vectors into 
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account, we arrive at the correct number of local operators. However an important 
question was left unanswered in [T7]: to identify precisely the descendants by BBS 
fermions and the usual Virasoro descendants. We shall address this issue in this 
paper. 

Being unable to solve this problem, the authors of [17] concentrated on the clas- 
sical limit showing that the description of the null vectors in terms of the fermions 
provides a new formulation of the classical hierarchy. Actually, they considered only 
right chiral descendants, so, the hierarchy in question was that of the Korteweg-de- 
Vries (KdV) equation. Another way of counting local operators was explained from 
a representation theory viewpoint in [TB] . 

Now we would like to discuss a seemingly completely different subject. Several 
years ago, the present authors together with H. Boos and Y. Takeyama started a 
joint work on correlation functions for the XXZ spin chain. We have been motivated 
by a strong feeling that the formulae known at the time were quite unsatisfactory. 
As a result, we found a fermionic description of the space of quasi-local operators. 

We believe our fermionic construction to be important, so, we would like to ex- 
plain its intuitive meaning. Generally speaking, our understanding of quantum field 
theory is limited because very few exact results are known about models with in- 
teraction. Our intuition relies too heavily on free fields, and this can be sometimes 
misleading. To give an example in the context of the XXZ spin chain, the free model 
is the XX spin chain arising at a particular value of the coupling constant where 
the model becomes equivalent to the lattice Dirac fermion. In this case, one diago- 
nalises the Hamiltonian introducing the creation-annihilation operators by Fourier 
transform. We parametrise the corresponding momentum as ip = log^^^. On the 
other hand, the same ^-dependent creation operators can be used for constructing 
quasi-local operators from a given one: we just take the (anti)-commutators and 
develop the result around the point ("^ = 1. So, the same construction with Fourier 
transform serves two different goals: diagonalising the Hamiltonian, and describing 
the space of quasi-local operators. It is important to understand that these two 
procedures are completely different for models with interaction. Moreover, before 
finding the fermionic description of the space of quasi-local fields we even did not 
know that the second procedure makes sense. Let us be more explicit about this 
point. 

The diagonilisation of the Hamiltonian for integrable models in general, and the 
XXZ spin chain in particular, is achieved by the Bethe Ansatz which is best under- 
stood in its algebraic formulation (ABA) [19]. The ABA can be viewed as a highly 
non-linear analogue of Fourier transform in the space of states. It does not explicitly 
introduce the creation-annihilation operators, and cannot be used for creating the 
local fields by adjoint action as it was possible for the free-fermion case. The real 
achievement of our works [201 EI] is that we were able to find a set of creation oper- 
ators b*((^), c*{() which act on the space of quasi- local operators creating this entire 
space from the "primary fields" by development around = 1. It may be said that, 
comparing to ABA, we have introduced another non-linear Fourier transform in the 
space of operators. In doing that we used the same algebraic structures as the one 
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used in ABA (quantum groups [221 ES] essentially) but in a more sophisticated way, 
including in particular the methods developed by 

Our fermionic operators have two notable features. First, they indeed act on the 
space of quasi-local operators, i.e., they respect locality. This was proved in the 
paper [21]. Second, the partition function of the the equivalent 6- vertex model, 
formulated on a cylinder with an insertion of a quasi-local operator created by the 
fermions, is expressed in terms of a single function uj{(,C,)- This was proved in the 
paper [2S]. We shall refer to the compact direction on the cylinder as the Matsubara 
direction. 

Our next goal was to take the scaling limit in order to arrive at the c < 1 CFT. 
After the transformation to fermions has been done this scaling limit is simple. In 
a certain sense we consider the fermionic construction for the lattice model as an 
existence theorem: the local operators are parameterised by parameters (, and the 
partition function with insertion is expressed in terms of u{(,C,)- To consider the 
scaling limit it suffices to describe it for u{(,C,)- Here the TBA-like equations for 
(^(CjO [26] are very useful. (We use the term TEA in a broad sense. Actually the 
techniques used in f2E[ [T] is that of ^7].) Comparing the scaling limit of u{(,^) 
with the CFT three-point functions, we find the relation between the description of 
local operators by the fermions and the one in terms of the Virasoro algebra. 

For the application to the sG model, we use an inhomogeneous XXZ spin chain, 
and obtain in the scaling limit fermions /32i-i) 72^-1) /^2j-i5 72i-i) which we call 
Boos-Jimbo-Miwa-Smirnov (BJMS) fermions. This construction is applied to solv- 
ing the longstanding problem of computing the one-point functions for the sG model 
on the plane [2] and on the cylinder [5]. To be precise, in addition to the BJMS 
fermions, certain fermionic screening operators are used to create the primary fields 
^Q+2mi^ their descendants from $q,. We shall not discuss them as they will 
be irrelevant for the purpose of this paper. Let us emphasise one more time that 
the relation between the BJMS fermionic descendants and the Virasoro descendants 
can be computed. This has been done in the quotient space by the action of the 
local integrals of motion up to level 6 in [T] and on the level 8 in [2H]. For studying 
the one-point functions it is sufficient to work in the quotient space. However, it is 
not quite sufficient for the goal of the present paper, we shall comment on this point 
soon. 

Though the BBS fermions and the BJMS fermions have been introduced by dif- 
ferent methods and for different reasons, there is a certain similarity between the 
two. This similarity motivated us to investigate the situation closely. On one hand 
we have the form factor formulae for the descendants written in terms of the BBS 
fermions. On the other hand we have local operators created by the BJMS fermions. 
Quite generally, if we insert such a local operator on the cylinder and take any eigen- 
vectors of the Matsubara transfer-matrix as asymptotic conditions, then the parti- 
tion function can be expressed in terms of a single function u{(,C,)- In particular, 
we can put an excited state to the left and the ground state to the right. Then the 
infinite volume limit in the Matsubara direction can be performed. It is clear that 
the result is nothing but a form factor. This provides us with the possibility for 
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comparison. To our great surprise, the BBS and the BJMS fermions are completely 
equivalent. In the multi-index notation (see flll.l9p ) the statement is this: 

(1.4) If 0^ = (3]+Wl+l*I-l*I-^a, 

(1.5) then L« = ^^d^-l+l^^l+l^-l+l^-'DV^^.V^t^Xj-Xl-MW , 

where |/| denotes the sum of the entries of the multi-index J. To be precise we have to 
add some fermionic screening operators in (11. 4p for m 7^ 0, but in the infinite volume 
they are irrelevant as explained in Section [Hi For us the precise identification (11. 4p - 
(II. 5p came absolutely unexpected. It demonstrates a remarkable self-consistency of 
the sG model: taking two complicated problems, that of computing the form factors 
of descendants and that of computing the one-point functions on the cylinder (at 
finite temperature) and going to the very bottom of them we find the same fermionic 
structure. 

The profit from the identification (II. 40 - 01. 5p is twofold. 

First, since the BJMS descendants can be quantitatively related to the Virasoro 
descendants, the form factors of the latter can be computed. There is one technical 
obstacle here: up to now we were able to identify the BJMS and the Virasoro 
descendants only modulo the action of the local integrals of motion. This was quite 
sufficient for one-point functions, but for form factors we would like to have the 
complete answer. The technical problem which one needs to solve for this goal is 
explained at the end of Section [3 

Second, in [17] the null-vectors are found in terms of the BBS fermions. So, 
using (ll.4p - (ll.5p we can identify them for the BJMS fermions. To keep the present 
paper within a reasonable size, we leave the detailed study of the null-vectors to a 
separate publication. Here we shall consider only the chiral null-vectors for the field 
^(2m-i)i^) m > 1 which corresponds in the CFT language to $i,2m. Let us give 
the simplest example which is the singular vector on level 2m. This singular vector 
is written as 

At the moment comparison can be made only modulo the local integrals of motion, 
but even with this simplification the Virasoro counterpart of (II. 6p looks really hor- 
rific. Using the formulae of [H |28] we find perfect agreement up to level 8. The 
fermionic formulae are also simple for other null vectors which are the descendants 
of the singular vector in the Virasoro language. We think that this nice simplicity 
is another evidence of the universality of the fermionic picture. 

The plan of the paper is as follows. In Section 2 we briefiy review the form factors 
axioms in application to the sG model. In Section 3 we discuss certain integrals 
which play a basic role for the description of the sG form factors. In Section 4 we 
present the formulae for the form factors. Section 5 introduces the BBS fermions. 
In Section 6 we briefiy discuss the BJMS fermions in the sG case. The origin of 
the fermionic description is the 6- vertex model as explained in Section 7. In Section 
8 we discuss the infinite volume limit in the Matsubara direction. In Section 9 
we present the main technical achievement of this paper, namely the computation 
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of the function u{(,^) for infinite volume limit in the Matsubara direction. The 
equivalence of the BBS and the BJMS fermions is established in Section 10. The 
null vectors for $i,2m in terms of the BBS fermions are discussed in Section 11. 
Section 12 is devoted to comparison of these null vectors in the fermionic and the 
Virasoro descriptions. 

2. Form factors axioms 

We are interested in form factors of the exponential fields $a(0) (11.21) and their 
descendants. In this paper we define form factors to be matrix elements of a local 
operator taken between excited states on the left and the vacuum state on the 
right, thus changing the convention of [7] where the opposite matrix elements were 
mostly studied. They are simply related by the crossing symmetry. We shall use the 
notation 0^ for descendants. These operators do not carry the topological charge, 
so, their form factors are non-trivial only in the case of an equal number of solitons 
and anti-solitons. 

We do not consider the breather form factors since they can be obtained as residues 
of soliton ones. So, we have 2n particles (n solitons and n anti-solitons) with ra- 
pidities ■ ■ ■ f32n- The form facotrs /©^(/Si, ■ ■ ■ , /32n) are meromorphic functions 
of these rapidities. They are vectors from the zero weight subspace of the space 
(C^)®^". The standard basis in the j-th tensor component is denoted as e^. 

The form factors are subject to three axioms. Formulating these axioms we follow 
the conventions of [7], namely, if two rapidities interchange we assume that the 
corresponding tensor components are permutated at the same time. 

Symmetry axiom. 

(2.1) %+i(/3,--/3,+i)/o„(/3i,--- - ,/32n) 

where Si^ is the soliton S-matrix [1]. Its explicit formula will be given later (110. 3p . 
Riemann- Hubert problem axiom. 

(2.2) foM, ■ ■ ■ , /32n-l, /32n + 2m) = e-^^'^^foMn. , P2n-l) . 

Residue axiom. 

(2.3) 27rZ reS^2„=/32„_,+^i/o, (A, ■ ■ ■ , /32n-2, hn-U P2n) = 

(l - e-^^"'^-^2n--l,l(/32n-l " A) " " " ^2n-l,2n-2 (/32n-l ' P2n-2)) 
X /Oc(/3l, ■ ■ ■ , /32n,-2) ® S2n-l,2n , 

where Sij = ef ^ ej + e~ ® e+. 

We change the standard basis of the tensor product to a new basis w^^'"' {(3i, - ■ ■ , (32n) 
described in [7], and express the form factors as 

(2.4) foMr-- ,/32n) = Z{(3u--- ,(32n) 
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where = {j | 1 < j < 2n, ej = ±}, and the sum over ej's is such that 
We have introduced an overall multiplier 

2n 

^(/3i,---,/32n) = ^nC(A-/3.)-e''"^' ^ 
The formula for the function ({(3) can be found in [7], 



2(1 - i^)7r2C(-7ri) 

We have set also 

/3/ = {Ai, ■ ■ ■ , A„}, if / = {zi, ■ ■ ■ ,^„,} . 

We use this notation only for symmetric functions, so, the ordering of the indices 
ii, ■ ■ ■ ,in is irrelevant. The main property of the new basis and of the function ({(3) 
is that 

C(A - - ■ ■ ■ , A, A+i, ■ ■ ■ ,/32n) 

= C(A+i - A)w^^^'-'^'+^'^"-'^^"(/3i, ■ ■ ■ , A+i, A, ■ ■ ■ , An) . 

Due to this property the first axiom is satisfied provided the essential part S^Ocni/^i- 1 A^ 
of the form factor is symmetric separately in /?/- and . It is well known that for 
this function one can write an integral representation. The integrals involved in this 
representation are quite remarkable, and we shall discuss them in detail in the next 
section. 

Before closing this section a remark is in order concerning the phase in the right 
hand side of (12. 2p . With solitons are associated some quasi-local fields, which inter- 
polate in- and out-states when time goes to =i=oo. The phase in (12.21) specifies the 
locality property of a given field 0^ with respect to these interpolating fields. The 
BBS fermions which will be discussed in section [5] act on a primary field and 
create fields which share the same locality property (i.e. the same phase) as for 



3. Integrals 

Consider the function 

2n 

x(o-| A, ■ ■ ■ , An) = Yl ^(^' l^i) ' 

i=i 

X{a, f3) = ie-^(^+'^-^ - /3 + f ) , 

where <^(cr) is defined in [7]. We shall not list explicitly the properties of v?(o"), since 
they can be read from those of the function x((t|/3i, ■ ■ ■ , An) which we are going to 
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give. In what follows we use the symbols 

s = e^", bj = e^^\ q = e"^i^, a = e"^^" . 

The function xi'^lPij ' ' ' i P2n) is a meromorphic function of a. For real /3j's it 
does not have singularities for > Im(cr) > — vr. It has the following asymptotic 
behaviour for a — )■ ±00: 

X(cr|/3i, ■ ■ ■ ,/32„) a:~(s|bi, ■ ■ ■ , b2„)X"(S'|5i, ■ ■ ■ , 52„) , 

where 

00 

(S|bi, ■ ■ ■ , b2n,) = 1 + ^X^(bi, ■ ■ ■ , b2n,)s"'' , 



X 



k=l 

00 



X+{S\Bu ■■■ , B2n) = 1 + 5^ X^{B^, ■ ■ ■ , S2„,)5-*^ , 

fc=l 

2n 00 

2:-(5|bi, • • ■ , b2„) = q" n (1 + ■ ■ ■ ' Ms') 

j=\ k=i 

2n 00 



j=l k=l 

The functions x^, (resp. x^, X^) are symmetric (resp. symmetric Laurent) 
polynomials of their arguments. They can be inductively computed from the func- 
tional equations 

(3.1) xicr + 2m)p{sq'') = x{(r)p{5q^) 

(3.2) xicT + '-^m)PiSQ) = xicr)P{-S) , 

where 

2n 2n 

(3.3) P{S) = l[{S-B,), p(5) = J](s - b,) . 

j=i i=i 

The normalisation of x(o"; /3) is defined by the identity 

1 



(3.4) x(c^, /3)x(f^, /5 + ^0 = ^ j,^(^2 _ 52) ' 

which is convenient for computing the residues fl2.3p . 
Consider now the integral 

(3.5) /a(/3i,--- ,/32n) = J xWu--- ,P2n)e^''da. 

R-iO 



10 M. JIMBO, T. MIWA AND F. SMIRNOV 

From the asymptotic behaviour of x(<^|/3i, ■ ' ' :/^2n), it is clear that this integral 
converges for 

2n 

(3.6) < Re(a) < — . 

In what follows we shall often omit the dependence of these functions on . . . , /32n, 
abbreviating them to la, xlc"); x'^is) and X'^^S). 

We want to continue the integral Jq, to the entire complex plane of a. To this end 
for any /c G Z let us introduce Laurent polynomials of the form 



which satisfy the identity 

(3.7) p(5q-') = m^'^) (s) + a-'p(s)n('=) (sq"^) - p(sq-2)n(^) (s) , 
where p(s) is defined in (13. 3p . For example 

(3.8) nW(5) ^ 



a-2q4n - 1 ' 

(2)r.^ - ^ / 1 ^ q'"(l-q') -2^^ ..... , ^ -1 



(3-9) n^^H^) - ^r,^ ^1 + ^;_.^-....-) ^-V.(b. . . . , b.„) . 

Here and after, a/(bi, ■ ■ ■ b2n) stands for the /-th elementary symmetric polynomial 
in bj's. Note that 



(3.10) res,2=,4„n(^)(s)— = -1, 

(3.11) resa2=„4n+4n''^^(s)^^ = —xtz^^. 

Using the functional equation (13. ip we can transform the integral defined in the 
original region (13. 6 p to 

(3.12) L = I x{a)e^''^^^da + J x{a)e^''n^'\,)da , 

R~iO r 



where the contour F contains all the poles of the integrand in the strip > Im(cr) > 
—271. It naturally splits into a sum of contours Tj as is shown on fig. 1. 
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Fig.l: Contours for analytic continuation of I a- 

The original contour T is split into a sum of the contours Tj . Each Tj surrounds a 
series of poles between the strip — tt > lm{a) > —2tx with the same real part as Pj. 

The equivalence with the original definition is established by transforming F to 
(-M - iO) U (M - 2ni - iO). But the right hand side of fl37T2D is analytical in 
< Re(a) < ^ + 2. So, we have managed to continue analytically to the strip 
^ < Re(a) < ^ + 2. From the formula dSS} we see that the result of the analytical 
continuation has simple poles at the points 2n + (2n + l)^-j^ with / > 0. We shall 
compute the residues at these poles later. 

Let us proceed. Another representation for /„ is possible for 2 < Re(a) < — + 2: 



(3.13) 



—' 



. m 



(2), 



R-iO 



p{sq' 



-da 



(2), 



s)da 



Indeed, the difference between (13.131) and (I3.12p is 

(3.14) I X(a)e^^-!^^i^i^|^^^ I x(a)e^^-(n(^)(.)-n(^)(.))rfa. 

M-iO r 

For 2 < Re(a) < — + 2 the integrals are well defined. Moreover, by deforming the 
contour as before, it is easy to show that (I3.14p is equal to zero. The right hand 
side of (I3.13P is well defined in 2 < Re(Q;) < ^ + 4. So, we have continued into 
the strip ^ + 2 < Re(a) < ^ + 4. The series of poles at 2n + (2n + /) ^ continues, 
and new poles at 2(n + 1) + 2{n + l)^-^ with / > appear. 

It is clear now how to continue further. It is equally clear how to continue it 
to non-positive Re(a). The final result is that la is a meromorphic function in C 
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with simple poles at 

a = 2(n + m) + (2n + /)i^, /,m>0; a = -2m-^l, l,m>0. 

Let us compute the residues. Consider /„ in the strip < Re(a) < 2 + 2n^. In 
this strip we have poles at the points a = 2n + (2n + 1)^-^, I > 0, because at these 
points m^^^ (s) and n*-^-* (s) have singularities. Obviously, 

Then, using this equality with f l3.10p in fl3.12p and transforming the contour F to 
-(-(-oo. A)) U (A - 2TTi, A) U (-oo - 27rz, A - 27ii) we find 

A 

res„=2n+(2n+o^^-^ = ' tL j x{^)e'''^'^^> da = -27rzX+ . 

A-27ri 

Next, consider Jq, in the strip 2 < Re(a) < 4 + 2nK Using f l3.13p . we compute the 
residues at the two series a = 2n + {2n + /) and a = 2n + 2 + {2n + /) of poles 
of n^'^^s) given by (13. 9p . For the former the result is the same because the term 
containing s disappears in the limit A — )■ oo. For the latter using (13. lip we obtain 

^6SQ=2n+2+(2n+/)i^-'a"^ ~ — /ITTZ^, . 

Continuing along the same lines we come to the following result. 

Proposition 3.1. The residues of 1^ at a = 2n + 2m+ {2n + 1)^-^ for 1,171 > are 
given by 

(3.15) res I^da = res ( S'X+{S)^^ res ( s'"x+(s) — 

a=2(n+m)+(2n+0i^ " S=oo \ S J s=oo \^ 5 

Similarly, the residues of la at a = —2m — l^-^ for /,m > are given by 

(3.16) res I^da = 1^ res f S~'X-(S)^] res ( s-'"x-(s) — 

a=-2m-l^ " S=0\ S J s=0 \ S 

Now we give the main definition. 

Definition 3.2. Consider two Laurent polynomials £{5) andL{S). We define their 
pairing {£, L)^ by the following two requirements: 

1. The pairing is bilinear. 

2. Ifi{s) =5^, L{S) = then 

(3.17) {i,L)a = Ia+2^_,i^i- 

The polynomials £ and L are interpreted as cycles and forms in the classical limit. 
Actually there are two possibilities. In the limit u ^ 1, i describes cycles and 
L describes one-forms. In the limit u ^ 0, i describes one-forms and L describes 
cycles. The situation is very much dual, and there is no preferred choice to call one a 
cycle and the other a form. So we will call i a q-deformed form and L a Q-deformed 
form. The next two propositions describe exact forms. 
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Proposition 3.3. For any Laurent polynomial z{s) we define 
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(3.18) Da[z\{5) = a-^p{s)z{5) - p{5q^)z{5cC) , 

and call it a (\- exact form. For any Laurent polynomial L{S) we have 

(3.19) {Da[zlL)^ = 0. 

Proof. It is sufficient to consider the case L{S) = SK Tlie procedure described above 
gives for any £(s) 

(3.20) (£,L)„= I Xia)e^^^^^-^^^da + J xi-)n{.)e^^^^^'^da , 

R-iO r 

where 

(3.21) p{sq-^)£{s) = m{s) + a-^p{s)n{sq-^) - p{sq-^)n{s) , 

and m(s) is chosen in such a way that the integral in (13.201) converges. This is 
possible to do with a comfortable margin. If we take i = Da[z] it is easy to see that 
(I3.2ip is satisfied with 

(3.22) m(s) = 0, n{s) = p(sq2)z(sq^) . 

Then the right hand side of f l3.20p vanishes since p(sq^) cancels the singularities of 
x((t) inside F. □ 

Proposition 3.4. For any Laurent polynomial Z{S) we define 

(3.23) Da[Z]{S) = Z{S)P{S) - AZ{SQ)P{-S) 

and call it a Q-exact form. For any Laurent polynomial £{5) we have 

(3.24) (£,D^[Z]), = 0. 

Proof. It would be sufficient to consider a simple case when there is no need for 
regularisation of the integral, and then continue analytically. However, in what 
follows it is more instructive to have a direct proof for any a. Figure 2 illustrates 
the proof. 
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71 



1-v 



-71"- 



-271 



1-v 



r+7i:i 

J V 



Fig. 2: Possible poles arising in the proof of Proposition [X^ 

Poles canceled by the factor P{S) are represented by ordinary crosses inside circles. 
The residues at the remaining poles depicted by boldface crosses cancel each other. 



From the definition we have 
(3.25) 



R-iO 



m[5j 
p{sq^' 



d(J+ / xi(^)e^''DA[Z]{S)n{s)da 



where m(s) and 77,(5) are polynomials satisfying fl3.2ip and are chosen by the re- 
quirement of convergence. Actually they may be different for different monomials 
in Z{S), but this does not matter for the following computation. 
Using the functional equation (13. 2p one easily finds that 



/ 



{i,DA[Z]), 



\ 



x{a)e^''P{S)Z{S) 



m[S) 
p{sq^' 



-dcr 



Y-iO M.+7Ti^~iO ) 

^ f\ 

J- J x{^)e^''P{S)Z{S)n{5)da. 

Poles arise from xi'^) ^"^^ We ask which poles are inside the difference 

of contours (M - iO) - (M + - iO), or T - (F + vri^^). The poles of the first 

integrand at a = (5j are canceled by P{S). So, the only remaining poles between 
M and M + vrz^^ are at a = f3j + 7ri(l — ^^[y^])' Similarly, the poles inside 
r — (r + i^^) are situated only at the points a = (3j — mil + [ir^])- C)n fig. 
2 the pole canceled by P{S) is in the circle, and the remaining poles are fat. It 
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follows from the functional equation (13.11) and the definition f l3.2ip that the residues 
at these remaining poles cancel each other. 

□ 

The last definition which we would like to give here concerns the multiple inte- 
grals. Consider k Laurent polynomials £i{s) , ■ ■ ■ ik{s) and k Laurent polynomials 
Li{S), ■ ■ ■ , Lk{S), and define the antisymmetric Laurent polynomials of k variables: 

£«(5i, ■■■5k) = {iiA- ■ ■Aik){5u ■ ■ -Sk), L^''\S,, ■■■ ,Sk) = (LiA- • ■ALfc)(5i, ■ ■ ■ , Sk) . 
We define 

(3.26) (£«,LW)„ = det((4L,)„),^^.^,...,. 

Then this pairing is generalised for all antisymmetric Laurent polynomials of k 
variables L^'^^ by bilinearity. 

4. Formulae for form factors 

Now we are ready to write down the formulae for form factors (12. 4p . Let us first 
prepare a couple of symbols. 

For a partition /" U /+ = {1, ■ ■ ■ , 2n} such that tl(-^~) = ttl-^^); define the poly- 
nomials 

Pi- (5) = n " = n ~ ' 

j&i- j&i+ 

so that we have p{s) = p/+(s)p/-(s). We set also 

where [ ]> signifies the polynomial part. Essentially following [7] we define 

(4.1) 4-U/+(5l' ^Sn) = {il-Ul+,0 A • ■ ■ A il-ul+,n^l) (Sl, " " " ,S„) , 

(4.2) ii-ui+A^) = ci~^{Pi~i5) (p7+,i(s) -p7+,i(sq2)) 
Formula (14. 2 p can also be rewritten as 

n-l 

(4.3) cj-ui+{t,s) := 5^(qH)"-Y,-u,+,,(5) 

1=0 

at . 2 . a"^q^t , , 

Now consider the ansatz (12.41) . As it has already been mentioned, the symmetry 
axiom (12. ip is satisfied automatically if 3^0cn(/3/- 1/3/+) is symmetric separately in 
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and in Furthermore, the Riemann-Hilbert problem axiom fl2.2p is satisfied 
if we set 



2n 



(4.4) 9a.„W,'\P,.) = e"- A" ' . (f^,, , , 

where L^"'^ = L^'^\Si, ■ ■ ■ , Sn\Bi, ■ ■ ■ , B2n} is an arbitrary Laurent polynomial which 
is anti-symmetric in Si^s and symmetric in -B/s. This statement is proved by a direct 
computation similar to the one given in [7]. 

While these two axioms concern form factors with a fixed number of particles, 
the third, the residue axiom, relates the Laurent polynomials L^"'^ with different n. 
Using the procedure of computing the residues familiar from [7] , one can reduce the 
residue axiom further to a simple system of recurrence relations given below (see 
(14.51) ). For later use let us formulate it in a slightly more general setting. 

Let c be an integer, and consider a sequence of Laurent polynomials 

lW = {L(''-)(5i,-- - - ,i?2„)}^n>o 

l — n — c 

which are anti-symmetric in S'j's and symmetric in -Bj's. We refer to L*^''"^ as the 
2n-particle component of L^*\ and c as the charge. 

Definition 4.1. We say that L^*^ is a tower of charge c if 

(4.5) L(''")(^i,--- - ,52.-2,5,-5) 



= {-lyB 1[{B' - S'^) ■ L('-i'"-i)(^i, ■ ■ ■ , S,„i|5i, ■ ■ ■ , 52„_2J 
p=i 

holds for all l,n > 1 with I — n = c. 

In the case c = 0, we also write L^*^'") as L^'^\ 

The statement is, the residue axiom is satisfied if the sequence {L'^")}^q in (14. 4p 
is a tower of charge 0. In other words, a tower of charge 0, or more precisely its 
residue class modulo Q-exact forms, represents a descendant field 0^. 

The most basic example of a tower of charge is given by the polynomials which 
do not depend on the parameters 5i, ■ ■ ■ , 52„: 

(4.6) M^")(5i, ■■■ ,Sn) = ■s^s''^■■■^ s^'^-' . 

This tower, denoted by Mq*\ is identified with the primary field whose vacuum 
expectation value ($q,) has been found in [T3] . 



5. BBS FERMIONS. 

We are now in a position to introduce fermions which create towers out of the 
tower Mq*^ for the primary field. In all formulas discussed below, we shall suppress 
the dependence on the parameters Bj. Since we deal with towers, however, their 
number 2n can vary from place to place. For that reason, in this section alone, we 
indicate the n-dependence by a suffix, e.g., Pn{S) = Yl^=iiS ~ Bj). 
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First we note that on the space of all towers there is an obvious action of the local 
integrals of motion. This is because the recursion relation f l4.5p is unaffected by the 
multiplication L^''") ^ /(J, J)L(''"), where /(/,/) is an arbitrary polynomial in 



(5.1) /2,-l,n = ^2,-l,n = J] ^ 



2n 2n 

(2i-i) 

k 

k=l k=l 



In the construction of towers we shall fully make use of the freedom ( 15. ip . We 
remark that the eigenvalues (15 .ip arise in the expansion 



(5.2) 
where 



oo ^ oo ^ 

(5.3) X„(Z) = ^-3T^"'''"'^2,-i,n , X„(Z) = Y ^-3T^''"'4-i." . 
Now, following [17], let us introduce the polynomial 

4z/ ei5i + e2S2 

ei,e2=± 

where the overall coefficient is introduced for future convenience. This polynomial 
is characterised by the following requirements: 

(5.5) Cn{Si, S2) is odd in 5*1 and even in S2, 

(5.6) degs,Cn{Su 52) = 2n - 1, degs,Cn{Su S2) = 2n - 2 , 

(5.7) Cn{Si, S2)\B2r,-i=B,B2-a=-B = {Si - B'^){SI - B'^)Cn~l{Si, S2) , 

with Ci{Si, S2) = —^{Bi + B2)BiB2Si. We note also that the total homogeneous 
degree of Cn{Si, S2) in Si's and B/s is 4n. 

Due to the definition through the pairing (14. 4p . adding Q-exact forms to compo- 
nents of a tower L*^''"-* does not change the form factors. A simple way of fixing this 
freedom is to restrict their degrees by 

(5.8) 0<deg5X(''")(^i,... ,50 <2n-l. 

If a is generic, any polynomial can be brought to this form by adding Q-exact forms. 

Let To,c denote the space of all towers of charge c satisfying (15. 8p . We introduce 
the action of fermions 

ro{Z) : V ^ ^0 ,c+i5 Xo{Z) : To,c — 7'o,c-i , 
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by defining 

(5.9) (V;*(Z)lW)^'+''"\^o,---,^z) 

= \ Skew5o,..,5,C„(^, S,)L^'^-\S,, ■■■,Si), 

(5.10) (xS(^)i^(*0^'^''"^(5i,--- ,5,_i) 

= jr^l {L^''''\Z,Su--- ,Si^,)-L^''-\-Z,S,,--- ,Si^^)) . 

Here g To,c, and 

Skew/(xi, ■ • ■ ,x/) = sgncr/(x^(i), ■ ■ ■ ,x^(/)) 

stands for skew symmetrisation. Because of (15.61) and (15.71) . ipQ^Z), Xo(^) preserve 
tlie degree condition (15. 8p and the recurrence relation (14. 5p . It is also easy to see 
that these operators mutually anticommute. 

Choosing representatives of restricted degrees (15. 8p is convenient to prove the 
completeness. Namely, following [IB] it should be possible to prove that ipQ^Z), 
Xoi^) create the complete set of solutions to the recurrence relations over the ring 
of hj-ii hj-i- 

Now we are confronted with the task of identifying these towers with local oper- 
ators. Notice that the coefficients entering the form factor axioms (I2.ip - (l2.3p are 
periodic in a of period 2(1 — z/)/z/. Hence for all integers m the axioms are the same 
for the fields 2^1^ ^-rid their descendants. One has to be able to identify the 
towers corresponding to all these fields. Consider the simplest case: the tower of 

2n 

(5.11) M^") = ($„+2mi^) ■ '5""+' A 52-+2 /\ . . . /\ ^2™+2m-l "Q p-m 

obviously corresponds to the form factors of the primary field 2^1^(0) • How 
does this fit into our description in terms of polynomials satisfying the restriction 
(15. Sp ? Obviously, for m > (resp. m < 0) the Laurent polynomial Mm^ contains 
degrees higher than 2n — 1 (resp. lower than 0). They have to be reduced using 
Q-exact forms. We need some general description of this procedure which looks 
completely random for the moment. 

The solution to the above problem was found in [T7j by considering the classical 
limit. In the classical case the variables Si, - ■ ■ , Sn turn into the separated variables. 
With our conventions the local classical observables are represented by polynomials 
L^Q^ of odd degrees in all Sj. So, the idea is to bring the towers created by iPq{Z), 
Xo{Z) to this form using Q-exact forms. We shall see that a very transparent 
structure will emerge. 

In order to make the fermions odd in Z, it suffices to multiply them by (15. 2p (that 
is, by making use of the integrals of motion). To make their action odd also in S, 
we modify C„(S'i, S2) by adding Q-exact forms in 5*2. Namely, we introduce formal 
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series C±,„(S'i, 6*2) such that 

(5.12) C-|-_„(S'i, 5*2) = C„(S'i, 5*2) (S^^ 00, mod Q-exact forms in S2) , 
and which are odd in both variables. This can be achieved by setting 

C±,n{Si, S2) = ^ ^ eie2-Pn(ei5'i)P„(e2S'2)r±(e2S'2/eiS'i) , 

ei,e2=± 

where 

00 

(5.13) r+(x) = - Y^i-xYi; cot f (a + ^) , 

1=0 
-1 

(5.14) T_{x)= ^(-x)'^cotf(a + i;). 

l=—oo 

Expand C±^n{Z,S) in Z^-^. Their coefficients of Z^-'"^ where j < ot j > n are 
Q-exact forms in S, and the remaining terms read as 



(5.15) C+,„(Z, S) = J2 ^^"~^^'^V2i-i(5) (mod Q-exact forms in S) , 

n+j 

P2j-liS) = ^ S'^''"V2j-l,2fc-l , P2j-l,2n+2j-l = ^ " COt |(a ' ^^""^ 



1/ 



(5.16) C_,„,(Z, 5) = Z^-' "'^j>2i-i('S') (mod Q-exact forms in S') 



P2,-l(5) = Yl ^'"~'''"%-l,2A:-l P2,-l,2n+2i-l = ^2n(5)' " ^ COt f (« - ^). 
k=l 

After these modifications we define new fermions as formal series at Z^^ = 00: 

(5.17) {^p'*{Z)L^*Y''''"\So,--- ,S,) 

= ^Skew5o,..,5,C+,„(Z, 5o)^(''"H^i, ■ ■ ■ , ^0 , 

(5.18) (^'*(Z)LW)^'+''"\5o,--- ,50 

^Skew5o,..,5,C_,„,(Z, 5o)i^^''"H^i, ■ ■ ■ , ^0 , 



VPn(^)Pn(-^) ^! 

where the right hand sides are understood as power series expansions in Z^^ . Define 
also X *{Z)y X by the same formula (15.101) as for Xo(-^)' replacing 

1 1 

Pni-z) ~^ ^fpj^)pjpz) 

and Taylor expanding in Z^^. 
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Finally we define ip*iZ), ip*{Z), x*{Z), x*{Z) by a Bogolubov transform 



(5.19) 
(5.20) 



i)*{Z) = e-ip'*{Z)e^-, ip*{Z) = eV*(^)e~", 
X*{Z) = e^x*{Z)e^^, x*{Z) = e^x'*{Z)e~^, 



where S acting on the 2?7,-particle component is 
dZ dX 



2mZ 27iiX 



C+,n{Z, X)^{Z)x{X) + C.,n{Z, X)ij{Z)x{X) , 



C±,n{Z,X) 



1 



^PniZ)P^i-Z)^P^{X)P^i-X) 



C±^n{Z, X) 



In the above we have introduced formally the annihilation operators ijj{Z), tp{Z), 
x{Z), x{Z), which are canonically conjugate to ip*{Z), 'ip*{Z), x*(Z), x*{Z), re- 
spectively and annihilate Mq*\ Obviously, the same operators are canonically con- 
jugated to ip*iZ), ip*{Z), x*{Z), x*{Z). That is why we did not put prime in the 
notation. Note, for example, that 



riZ)=i;'*iZ) 



dX 
2mX 



C+,„(Z,X)x(X). 



The modification fl5.19p . fl5.20p is useful for separating the role of the two chiralities, 
as it will be explained shortly. 

With the definition of ip* (Z) , ip* (Z) , x*{Z),X*iZ) given above, let us write down 
their action on the tower Mq*^ . Suppose k — k' = I — n. 

(5.21) (riZi) ■ ■ ■ r{z,)r{Zp+i) ■ ■ ■ riZk) 

X X*{Xk') ■ ■ ■ x*{X,^i)x*{X,) ■ ■ ■ x*{Xr)MtY'"\Su ■ ■ ■ , 5,) 



\kk' 



n ^Pn{Z,)Pr,{-Z,) fi VPn{X,)P,{-X, 

j=i i=i 



A S 

e D 



where A, T), G and T) are respectively kxk',kxl,nxk' and n x / matrices: 
/ ••• C+(Zi,X,+i) ■■• C+(Zi,XfcO\ 



A 











C-{Zp^i,Xi) ■■■ C_(Zp+l,Xg 



C+{Zp, Xg + l] 







C+{Zp, Xy 





V c_(Zfc,Xi) 



C-{Zk-, Xq) 







/ 
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C+{Zp, Si) 
C~{Zp^i, Si] 



\ C-{Zk, Si] 
Xi ■ ■ ■ Xk' 



C+{Zi,Si) \ 

C+{Zp, Si) 
C-{Zp^i, Si) 



C4Zk,Si) J 
(Si ■■■ Si 



2n-l 



2n-l 



2n-l 



Si 



2n-l 



The Fourier modes are introduced by 

oo oo 

(5.22) r{Z) = 5^Z-2^-+V2Vi' X*{X) = Y,X~'^^\l^_, 



(5.23) 



r{z) 



oo 



X*{X) 



oo 



A similar determinant formula is obtained for ipl^Zi) ■ ■ ■ ipQ{Zk)xo{Xk') ■ ■ ■ Xo{Xi)M^\ 
for which the yl-part is 0. For completeness we give the Fourier decomposition for 
the annihilation operators 



^{Z) = Y,Z''-'^2,-i 



oo 



2^2j+l 



X2j-1, 



X{Z) = J2Z''^'X2,-1, x{z) 
and the non-vanishing commutation relations 

[V^2j-l,^2fc_l]+ = [X2i-l,X2fc-l]+ = [^2i-l,^2fc-l]+ = [X2j-l,X2fc-l]+ = ^j.fc • 

Let us explain the formula f l5.2ip using physicist's terminology, choosing k = k' 
for simplicity. Suppose we create a new tower starting from the primary one by 
application of several operators 4'2j-iy X*2j-i^ i^2j-i^ X2j-i- Let us take n ^ j for 
all indices j of the operators involved. Then the polynomial Mq^\Si, ■ ■ ■ , S'„) = 
{^a)S A S'^ A ■ ■ ■ A S*^""^ can be considered as Fermi zone with two ends situated at 
S and S*^""^. Roughly speaking the operator X2j-i creates a hole at 5'2»^-2i+i^ X2j-i 
a hole at S'^^~^, '^2j-i ^ particle at 5'2"+2i-i^ ^^^^ '^2j~i ^ particle at S~'^^~^^. 

The tower for the primary field is represented by the pure wedge product 04.61) . 
Those for the descendants are linear combinations of pure wedge products with 
symmetric Laurent polynomials in Bj as coefficients. We discuss which pure wedge 
products appear in the linear combination. Suppose aj, bj G 2Z>o — 1, and we pick 
up 

(5.24) (^:^ ■ ■ ■ ■ ■ ■ raM. ■ ■ ■ xi,X,, ■ ■ ■ xlM*^) ^"'"Vi, ■■■,Sn) 
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from dSJH) taking the coefficient of n'=i ^7"^ Tl.tp+i ^? n'=i ^/'^ n'=,+i - 
We assume that ai > . . . > Op, a^+i < . . . < ak, bi < . . . < bg and bg+i > . . . > b^. 
Moreover we assume that n is large enough so that 2n — bg > b^+i, too. Note that 
m = p — g is the excess of particles at the right end of the Fermi zone, and — m is 
the excess of particles at the left end. We call m a weight. A basic example of a 
tower of weight m is (15. lip . 

The procedure of taking the coefficient is as follows. First we take the coefficient 
of Zj (1 < J < p) and Zj^ (p + 1 < J < A;) in the ffist k rows. Second we 

take the coefficient of Xj {1 < j < q) and Xj^ + 1 < j < /c) in the ffist 

k columns. For example, take 1 < j < p and consider the coefficient of Zj ""^ in 



C+^n{Zj ,X)/ ^Pn{Zj)P,,{-Zj ) . We expand this for Z. — 7- oo. The leading term is 
obtained from the approximation 



the coefficient being paj{X). The sub-leading terms are multiples of p2i-i{X) where 
21 — 1 < ttj. Therefore, if we consider a partial-ordering of pure wedge products by 
the natural ordering of aj, the sub-leading terms has coefficients of lower order. It 
is the same for j> + 1 < j < A; and the column expansions. Collecting the coefficients 
in the leading order we obtain a matrix, whose {A, S) part is given by 

(5.25) 





Q Pai,bq+i 


■ Pai,bk PaASl) ■ 


■• PaASn) 


\ 




Pap, 6,-1-1 


■ Pap,bk Pap{Sl) ■ 


■■ Pap{Sn) 




Pap+i,bi 


Pap+i,bq ^ 


^ Pap+i{Sl) ■ 


■ ■ Pap+i (Sn) 




\ Paf^M 


Pak,bq 


Pa,iSi) ■ 


■■ Pa^Sn) 


/ 



Consider q + 1 < j < k and set bj = 2/ — 1. The l-th row in the (C, D) part gives 
rise to 



,1-st <j-th (g + l)-st j-th fc-th u. u 

(5.26) (0 ■■■ ... 1 ■■■ Q ■■■ S'n 



Similarly, for 1 < j < g, set 2n — bj = 21 — 1. The l-th row in the (C, 2)) part gives 
rise to 

(5^27) ('o' ... T ... 'o" "T* ... 'osr-- ■■■ Si"--') 

If 2/ — 1 ^ . . . ,bq,2n — . . . ,2n — bk} the l-th row in the (C, D) part gives 
rise to 

l-st g-th {g + l)-st fe-th , „, , 

( ■ • ■ ■ ■ ■ ^1 • • ■ sl^-^) 
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Subtracting certain multiples of rows f l5.26p . (15.27^ from the part f l5.25p . we 

change the latter to 



(5.28) 



(0 (p.,(5.)-EUll>a.,.,^^0.-,.. 

fc=i . 

(^Pa,{Sk) -T.LlPa.MS, 



\ 



2n~bi 
k 



fe = l,...,n / 



In physicist's terminology fl5.26p implies the existence of a hole at S^j , and fl5.27p 
hole at S'^^~^K On the other hand we obtain a row from (15.280 by expanding it in 
the same power in Sk'- 



l-st 
Pa,,2l-l{, 



fc-th 



Si 



2l~l 



n ) 



;i<j<p; 2/-l^{6,+i,...,6„}) 



or 



l-st 
Pa,,2l~l{, 



A;-th 



Si 



2n-2l+l 



^2n-2;+l) (p + 1 < J < fc; 2/ - 1 ^ {6i, . . . , b,}) 



The former creates a particle at ^ if 2/ — 1 > 2n, and fills the hole at 5*^" if 
2/ — 1 = 2n — hj for some 1 < J < Q'- The latter creates a particle at 5"^"^^'^^ if 
2/ - 1 > 2n, and fills the hole at S^^ if 2n - 2/ + 1 = hj for some q+l < j <k. It is 
important that the change caused by filling a hole by a particle occurs only if both 
particle and hole are near the right end, or both near the left end. Therefore, the 
weight m is invariant by this change. In other words, one can define the weight m 
of a tower without ambiguity. This is the effect of the Bogolubov transform (15.190 . 
(I5.20p which introduces the non-trivial block A. 

From (I5.15p . (15.161) it is easy to see that, up to a sign, the leading term in (I5.24p 
is given by 



2n p k 

(5.29) ($.)(^)^'n^rn^°tf(«+v) n ^oil{a-^-^)■s-'^^^■■■^s~'^^-^ 

A 5 A ■ ■ ■ (5^"=)^ ■ ■ • [s'-'+'y ■ ■ ■ (^2"-''^)^ . . . {s^^-'^y ■ ■ • A 5'"-^ 

^2n+ap . . . ^2n+ai 



Here {S^Y represents a hole. 

Our goal is to identify this fermionic picture with the content of local fields in 
the sG model. Here we do it qualitatively, leaving the quantitative identification to 
later sections. 
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Returning to the primary fields we find 



(5.30) ri ■ ■ ■ r2m-ixim-i ■ ■ ■ xi 



/ffS \ ™ 



(5.31) ^r-'-V^L-iXk-i-'-xI 



02m- 

^^"^ , ■ i-r TT cot f (ai^ - (2j - 1)) M 

V J — i 



m ■ 



In general, it is now clear that the space of descendants of Mq*'' by fermions, together 
with the action of the local integrals of motion, has the same character as that of 
the space of fields in CFT: 



V2i^m®y 



m=—oo 



where V„ denotes the Virasoro Verma module with central charge c = 1 — and 

2 

highest weight A^^ = 



:a(a — 2). 



We remark that often it is convenient to change the definition slightly and work 
with fermions given in (I5.2ip wherein the block S is replaced by 

/C{Z,,Si) ■■■ C{Zi,Si) 
(5.32) S ^ : : 

\C{Zk,S^) ■■■ C{Zk,Si) 

The resulting fermions have mode expansions fl5.22p . (15.231) in odd degrees. Acting 
on the primary field, the Fourier components generate towers which coincide with 
( 15.2ip modulo Q-exact forms, but satisfying the restricted degree condition ( 15. Sp . It 
should be noted, however, that we cannot replace C± to C in the block A because 
their second arguments are not the integration variables. 

6. BJMS FERMIONIC DESCRIPTION OF SG MODEL 

In the paper [1] we gave a fermionic description of CFT. This description has an 
advantage of being compatible with the integrable perturbation. As we explained in 
[21 [3] the logic of Perturbed Conformal Field Theory (PCFT) [H |29] implies that in 
a generic situation (z/ and a irrational) there is a one-to-one correspondence between 
the local fields before and after perturbation, i.e. they consist of primary exponential 
fields and their descendants created by two copies of the Virasoro algebra with 
generators 1^, 1^. Following [T] we consider the case 

< «< 2. 

Equivalence to the fermionic description implies that the local fields are created 
from the primary field by the local integrals of motion i2fc-i, hk-i and the fermions 
P'*2k-ii 72fc-i; /^2fc-i; 0'2fc-i- is explained in [HEH] how to recalculate the result of 



FERMIONIC STRUCTURE IN SINE-GORDON MODEL 



25 



acting with the fermions in terms of the usual Virasoro descendants. Unfortunately, 
the techniques of [H [2E] allow us to do that only modulo the action of the local 
integrals of motion. At the end of the next section we give some explanation to this 
point. This was not a problem for the papers [2], [1] where the one-point functions 
in the sG model were calculated, because they vanish on the descendants created 
by the local integrals of motion. But in the context of the present paper this is an 
unpleasant restriction. We hope that the technical difficulties behind this problem 
will be resolved in the future. The origin of these difficulties will be explained in 
the next section. So, for the moment we work with the space 

a ^ a 5 



oo 



where 

oo 

fc=l k=l 

Let us emphasise one more time that the operators 12^-1, 12^-1, (3*2k-ii l*2k-i^ ^Ik-i^ 
72fc_i are defined on Va and the restriction to V^"° is due to purely 
technical reasons. 

To be precise the space (g) allows a fermionic basis: 

(6.1) /3^^j+7j-7/-<^>a(0), 
where #(/+) = #(/"), #(/+) = #(/^). We set generally 

Pi Pai ' ' ' Pap 1 Tap ' ' ' O'ai i Pi Pai ' ' ' Pap i ^ I O'ap ' ' ' Tai ' 

for / = {ai, ■ ■ ■ , ap}, a\ < a2 < ■ ■ ■ < a-p ■ 

On the other hand the space can be realised as result of acting on 

the primary field $0 by even generators of two chiral Virasoro algebras l_2fc, l-2fc- 
The formulae relating the fermionic basis with the Virasoro basis can be found up 
to level 8 in [H |28] . The fermionic basis was used in [H [3] because it allows one to 
compute the one-point functions. 

We use the same convention for multi- indices as in [H |3] : 

if J = {ji, ■ ■ ■ , jp} then aJ + b = {aji + 6, ■ ■ ■ ,ajp + b} , 
I{m) = {1, 2, ■ ■ ■ , m}, Jodd(m) = 2J(m) - 1 . 

The operators (3i*2k^i, l*2k-i ^"^^ /^2fc-i) l*2k-i combined into the generating 
functions 

00 00 

(6.2) /3*(A) = 5^A-^/3;,_„ 7*(A) = J]A-^7;.-i, 

fc=l k=l 

00 00 

^*(A) = x'^Kk-i. 7*(A) = ^"^7^.-1 • 

k=l k=l 

In a weak sense these series describe the asymptotics of holomorphic functions at 
A — !■ 00 and A — ?■ for right and left chiralities respectively. When expanded at the 
opposite points, the same holomorphic functions create another interesting set of 
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operators (actually a slight modification is needed [HIS]). This allows us to define 
the operators 



oo 

(6.3) 

^screen Tscrccn,fc' /^screen 

(A) = J2 A"-''^ screen,/; 

k=l k=l 

We define as usual 
(6.4) 

Tscreen,/(m) Tscreeriim Tscreen,! ' /^screen, 7(m) /^screen,! /^screen, m 

Then the m-fold screened primary field is by definition 

m 

(6.5) <l>(-)(0) =z-M''"n^°*f (2j-a) x^:_^,(^) 7™,/M'^'a(0) , 

where the multiplier in the right hand side is introduced for convenience. 

^quo 
a+2m 



It is explained in [3] that 'V''"° ^ V^"°2^ 1-^ for m > can be constructed as 



(6.6) /3^^/+7/-77-$.+2,.i^(0) 

-C^(a)/3^+2^^;-+_2^7/-+2™7l-2™/3L.M7L.M^i'"H0), 
where for negative indices we set 

(6.7) 7-a = -^a(«)/3a, = "4(2 - a)7a , ta{a) = cot ^{au + a) 
the only non-trivial anticommutation relations are 

and 

/3,$(r)(0)=0, 7,$M(0) = 0. 

In particular, 

(6.8) $a+2™i^(0) - C^(«)/3L.M7L.M'^i"'H0) , 
where (see [3]) 



(6.9) C^(a) = ^^^(-^i.)- J] tan |;(az/ + 2j - 1) . 

^ i=i 

The fermionic screening operators /9screen.2i-i) 7scrcen,2j-i P^^J important role in 
the CFT computations. However, in the sG case when we consider the one-point 
function on the cylinder [3] their contribution completely factors out. This results 
in the possibility of the identification in the weak sense: 

(6.10) $M(o)^$^(0). 

This is a general fact for the sG model. We shall comment on it later. 

Considering the formulae (16. Sp . (16. 9p . fl6.10p we find an amazing similarity to the 
formula (15.111) . This similarity was one of starting points for our present research. 
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Recall that up to now we considered only the case < a < 2. According to (16. 6p 
it is actually sufficient to consider the fundamental domain 

(6.11) 0<a<2^. 

We excluded the point a = from the fundamental domain for the following 

reason. An important role in the definition of the BJMS fermions is played by the 
function A'^^ipQ{(,a) discussed in the next section and its asymptotics (I7.18p . For 
a = 2^^ the asymptotic formula at C — )■ breaks down, and a logarithmic term 
appears from the j = 2 term in the ffist sum and the j = 1 term in the second sum. 
So, the case a = 2^-^ requires a special treatment. 

7. SiX-VERTEX MODEL AS ORIGIN OF BJMS FERMIONIC DESCRIPTION 

We introduce fermions for c < 1 CFT and for the sG model taking the scaling 
limit of the homogeneous and inhomogeneous six- vertex model on the cylinder. This 
is explained in details in the paper [3], so, we shall be very brief here. But we would 
like to repeat some facts from the paper [25] which is the cornerstone for all of our 
recent researches. 

Let us consider the inhomogeneous six-vertex model on the cylinder. We count 
the sites in the space direction by indices j, and in the Matsubara direction by m. 
The model can be inhomogeneous in both directions with the parameters of inho- 
mogeneity (j and Tm. The Boltzmann weights are combined into the L-operator 
Lj^miCj/'^m)- Using the equivalent XXZ spin chain language we introduce the op- 
erator g2aS{o)Q ^j^]-^ Q being local and S{k) = \ Yl'j=-oo'^j ■ Then we consider the 
partition function as in the following figure: 



Space 



<N a+K| 

\ — , I 




Fig. 3: Partition function of the six-vertex model on a cylinder. 

To each site of the lattice is attached a local Boltzmann weight Lj ^n- On one horizon- 

3 

tal line a background field q'^'^ is applied. On the same line, one allows in addition 
a local dislocation in the middle accompanied by an extra field q'^'^'^ extending to the 
left, representing the insertion of a quasi-local operator q^°'^^^^ . 

The number of sites in the Matsubara direction is denoted by n. The boundary 
conditions in the space direction are given by the eigenvectors {N_,k, + a\ and 
\N^,K,) of the Matsubara transfer-matrices with twists k -\- a and n, and N± are 
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the labels for different eigenvectors. In this formulation the number of sites in 
the space direction is not important as far as the operator fits into the picture 
fig. 3. The same computation can be used in a different situation. Put arbitrary 
boundary conditions on the right end of the cylinder, and allow it to grow to the 
right infinitely. Then the partition function will look for the eigenvector \0,k) with 
the largest eigenvalue of the Matsubara transfer-matrix (ground state). Certainly 
one should not be extremely unlucky which means that the vector representing the 
boundary conditions should not be orthogonal to the ground state. 
The partition function fig. 3 is written as 

(iV_, K + a|Trs (Ts.Mg^'^+'^^^^^O) \N+, «:) , 

where S stands for space, M stands for Matsubara, 5* is the total spin in the space 
direction, Ts,m is the rectangular monodromy matrix in the tensor product of the 
space and the Matsubara Hilbert spaces. 

Let us make a digression. We discuss a connection to form factors. Consider two 
quasi-local fields with opposite twists, q~'^°'^^^^Oi and q'^"^^'^^02. For simplicity we 
take K = and denote corresponding ground state by |vac). In Section 8 we will 
see that in the infinite volume limit in the Matsubara direction, the K-dependence 
is dropped. In order to compute the correlation function 

vac) 

we have to glue two partition functions of this kind and sum over the intermediate 
states as in fig.4. 




vac> 



Fig.4: Form factor decomposition of the two-point function on the lattice. 



This is nothing but the form factor decomposition in the Matsubara direction. We 
will come back to this subject in Section |9l 
In the paper [25] we computed the functional 

2K5+2aS'(0)( 



{N^,K + a\ 


Trs (T5,Mg''^^+'"^(°)0) 






|Trs (Ts,Mg2«5+2aS(0))|JV^^^) 



(7.1) ^«{gi^6+2a6(U)Q| 

The dependence of on the particular eigenvectors of the Matsubara transfer- 
matrices is not explicitly exhibited. Obviously (17. ip is not a form factor but rather 
a ratio of the form factor of a descendant field to the form factor of the primary 
field. This should be remembered. 
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The key objects to the computation of were introduced in the paper [21]. In 
this paper it was shown that the quasi- local operators g^^-^C^)© can be created from 
the primary field by action of one bosonic and two fermionic operators: t*((^), b*((^), 
c*(C). More precisely in the homogeneous case the quasi-local operators g2a5'(o)Q ^^^^ 
created by t*, b*, c* defined through 

(7.2) 

oo oo oo 

t*{c) = Y.(^'-irx, h*{c) = j2{e-ir'h;, c*(c) = 5^(c'-i)^-^c;. 

p=i p=i p=i 

In particular t^/2 is the right shift by one site along the lattice: 

(7.3) itt(g2"S{0)) _^2aS(l)_ 



The completeness of this basis is proved in [30 
The main theorem |25| states that 



(7.4) Z^{t*{r^,) ■ ■ ■t*(r/„Ob*(Ci) ■ ■ ■h*{Ck)c*{^k) ■ ■ ■ c*{^^){q'^'^'^)} 

m 

= n2p(^.)-det(a;(C,0))i,,=i,..fc, 



where 



Tk(Q 

and Tl(C) and Tr(C) are respectively the eigenvalues of the left (with twist k+q) and 
the right (with twist k) Matsubara transfer-matrices on the eigenvectors {N_, K + a\ 
and |A^_|_, n). The function a;((^, ^) was defined in [25] through a number of conditions 
as a quantum deformation of the normalised second kind differential on a hyper- 
elliptic Riemann surface. For computational purposes it is useful to consider an 
equivalent, alternative formula [2Sj which is explained below. 

Consider the eigenvector |iV_|_, k). The corresponding eigenvalue is defined by the 
function a{Q which satisfies the DDV equation 



(7.5) 



log a(C) = + log - I K{C/0 log (1 + a(0) ^ 



where the cycle 7 goes around the zeros of the eigenvalue of the Baxter operator 
Qr{C) (Bethe roots) in the clockwise direction, as opposed to all other contours. 
The functions a{() and d{() are defined by using inhomogeneous parameters in the 
Matsubara direction. The formulae for them can be found in [25]. We would like to 
emphasise that a(.^) depends only on the eigenvalue of the right transfer-matrix. 
For completeness let us write the Baxter equation, and the definition of a(C) 

TKiOQRiC) = diOQRiCq) + aiOQRiCq-') , o(C) - ^^^^^^^^^^ 



aiOQKiCq-') 
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Recall also that CQrIC) is a polynomial in The kernel in the integral equation 
is defined through a more general one 

K{O = K{C,0), K{C,a) = ^AcMC,c^), 
where we use the symbols 

Ac/(C) = fiCq) - fiCq-') , ^o(C, a) = C 



In the papers [251 [1] we used instead oitpQ^C,, a) the function ^/;(^, a) = ■ipo{C, a) — 
As it has been shown in [26] on the finite lattice, either of these two functions can 
be used in the equations defining without changing the final result. So, we 

can choose the one which is more appropriate for the scaling limit. 
In order to define oj{CyO have to introduce the operation S'^ 

s^fiO = fiCq) - p(C)/(C) , 

and the convolution 



with measure 



(l + a(,))- 

The fundamental role is played by the dressed resolvent 

(7.T) -RfJress -Rdress ^ ; 

where stands for the integral operator with the kernel K[C,/S,,ot)- Introducing 
two more kernels 

(7.8) /Mt(c,o = h 5c^o(c/e,«), /nght(c,o = 55>o(c/e,«) , 

we define (^(C^O by 

(7.9) = (/left ^ /right + /left * -Rdress * /right) (C; — Wo(C; 

where 

(7.10) a;o(C, = -5^ Vo(C/e, «) • 

The definition of the "primitive function" is ambiguous. Using Z = C^/", H = 
77^/*^, X = we shall define once and forever [1] 



oo 

(7.11) Vo(C/e, «) = ^VP J MC/V, «) 



H drf 
H + X 27rzr/2 
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where the principal value is taken with regards to the pole at rf = By this 
definition, we have 

oo 

(7.12) A-Vo(g^'C/e, «) = ±^77^ + ^ [ MCe^'Vv, «) ^ "^"^^ 





Here the left hand side is the analytic continuation of A'^^%jjo((, a) in the variable 
^ log C from Imi log C = to Im^ log(g^^C) = ivrz. 

A physically interesting situation occurs when we take the scaling limit which 
implies in particular n — )■ oo. Here we consider two cases. We shall be very brief 
because detailed explanations are given in [3]. 

1. Chiral CFT. 

In the homogeneous case Q = 1, = I, the scaling limit 

(7.13) n^oo, a^O, na = 27rR, ( = X{CaY , 

{R and A are finite) describes the chiral conformal field theory [1]. The important 
constant C is given by 

C= y^/, , T{i,)^ . 

2 (^) 

2. Sine-Gordon model. 

In the inhomogeneous case 

(7.14) Q = Ct''\ r^ = Cr'^ 







the scaling limit 

(7.15) n oo, a -> 0, na = 27iR, (q oo, = {C ay , 

[R and i-i are finite) we describes the sG model on a cylinder of radius R and the 
coupling constant n (see (11. ip ).. 

Let us briefiy recall the scaling limit of the operators t*((^), h*{() and c*{() in 
the inhomogeneous case. The construction is not trivial, so the interested reader is 
referred to the papers [21 [3] for details. The local operators in the inhomogeneous 
case are created by coefficients of b*(^), c*(^) developed into seria around = 
Co and = Co"^- From these seria, after certain Bogolubov transformation, we 
obtain the operators b+*((^), c"''*(C), b~*(C), c~*{(). Their expectation values are 
the determinants fl7.4p wherein the "two-point correlators" given by 

(7.16) Z„{b+*(C)c+*(O(g'"^(°0} = a;(C,0, 
Z„{b+*(C)c-*(O(g^"^(°0} = uiCO + cooiCO , 
Z„{b"*(C)c+*(0(g'"^(°))} = a;(C,0 + ^o(C,0 , 
Z„{b-*(C)c-*(e)(g2"^(°))} = u;(C,e). 
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The scaling limit is dictated by the behaviour of the functions p, u and ujq{C,,^). It 
can be shown (see Appendix A of [T] for very similar conclusions) that in a very 
general setting 



(7.17) lim logp(C) ^ Vpe(2j-i)C" 

scaling log C^^eoo ~^ 



scaling logC^eoo 

log^— >e'oo 

where e, e' = ±. It remains to consider uo{(,^). Recall the definition fl7.10p . Since 
we know the asymptotics of p(C) the only non-trivial thing is the asymptotics of 
A'^'^ipo((, a). For studying the asymptotics it is convenient to rewrite this function 
in the form of a Mellin transform: 



oo 



(7,18) A,-V.„(C.a)4/r + f)-i^^^^;(^* 

— OO 

oo oo .. 

c^oo '^'^ ^ ^'^^ sm7rz/(a — 2?) 

oo oo 

~ f V(-1)^C^ cot f (az/ - j) + i V ^. 

C^o^^^^ ^ ^""^ ' 2^ sin7rz/a + 2? 

i=i i=o ^ 

Motivated by the above formulae we have conjectured [2], [3] that the following 
scaling limit exists for the operators b^*((^), c^*((^): 

^ /3*(K) + (i + 0({i4))^™(C//^), 

s>oo 

^ y(A^C) + (i + o({i4))7™(C/M), 

s>oo 

^ ^*(c//x) + (i + o({i.}))/3:_(K), 

^ r (C/a*) + (1 + 0({i4))7™(K) , 

where 0({i^,}), 0({i*}) stand for descendants created by the local integrals of mo- 
tion. We have to apologise for forgetting to write these terms in |3] . There are two 
explanations for that. First, in [3] we considered the one-point functions for which 
the descendants by the local integrals of motion are irrelevant. Second, most im- 
portant, these terms are absent in the asymptotics for the chiral CFT which serves 
to normalise the fermionic operators. Their appearance is an artefact of mixing two 
chiralities, and we never use them. 

The operators in the right hand side are the unique [21 |3] operators in the sG 
model which provide deformation of the CFT operators (16. 2p . (16.31) . The prescrip- 
tion for powers of the spectral parameter are the same: they are dictated by the 
scaling dimensions. Looking at the asymptotics (I7.17p . (I7.18P one concludes that 
the pairings of the screening operators come only from a;o(C;0- Moreover, these 



ib+*(C) 


^ /3^*(C) 

scaling ' \ J / 


ic+*(C) 


— > 7+*(C) 

scaling ' \ j / 


ib-*(C) 


scaling ^ ( } 




— y '7~*(C) 

scaling ' \j/ 
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pairings are diagonal: /^scrccn.j couples only with 7*creenj 7screen,i couples only 
with /^screenj- The Corresponding pairings are easy to compute, and as a result we 
come to the conclusion that the screening operators can be ignored provided we 
postulate the formula (16.101) . 

An important generalisation of the above construction was proposed in [T1I3]. We 
have three parameters: a, the twist K + a for the left Matsubara transfer- matrix, and 
K for the right Matsubara transfer-matrix. By introducing the screening operators 
in the lattice construction it was shown that in the scaling limit we can achieve the 
emancipation oi n + a, namely, it can be replaced by an arbitrary parameter k'. The 
basic formula fl7.4p and the definitions of p{() and u{(,^) remain valid in this case. 
A particularly nice situation occurs for n' = n which implies p{Q = 1. This case is 
physically relevant because it describes the one-point functions on the cylinder for 
the sG model [3]. But also this case is simpler from the technical point of view, and 
it allowed the quantitative investigation of the equation for in the CFT case 

[H [2H]- The price to pay for this simplification is that we can identify the action 
of descendants created by fermions with those created by the Virasoro generators 
only modulo the action of the local integrals of motion. It would be important to 
remove this technical obstacle. In Section [9] we shall see the first example of an 
exact solution for at p{Q ^ 1. 

8. Infinite volume in Matsubara direction. 

Consider the limit /2 — ?■ oo which corresponds to the sG model on the plane. In 
this case the scaling limit and the limit of the infinite volume for the lattice model 
commute. In other words the correct result for the functions p and u will be achieved 
if we consider the limit n — )■ oo directly in the six-vertex case without the rescaling 
(I7.13P or (17.151) . This is so well known that we have a problem with making proper 
references. However, there is one less known point which we would like to underline 
in this section. 

The results of this section is used for both left and right transfer-matrices, so, we 
shall denote the twist by ^, assuming that it is real. In the limit n — >■ oo, the zeros 
of Q(C) densely fill M_|_ in the plane of . Obviously the solutions to the equation 

(8.1) a(C) + l = 0, 

are the zeros either of (5(C) or of T(C). For the ground state all the real positive 
solutions to the equation (18. ip are the zeros of (5(C)- Fo^' the excited state a finite 
number of real positive solutions are the zeros of T(C). They are conventionally 
called holes, and correspond to solitons in the sG language. In the sector of zero 
spin the number of holes is even (this is our In). It is necessary to introduce 
the same number of complex zeros of (5(C) order to compensate the spin. The 
positions of these complex zeros are defined by the positions of the holes through 
the Higher Level Bethe Ansatz (HLBA) equations [31]. The HLBA equations have 
(^") solutions which count different isotopic structures for solitons-antisolitons. One 
more property of these complex solutions is that they do not contribute to the local 
integrals of motion taking care only of the isotopic structure. We shall discuss this 
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in more details in Section [TOl There may be an additional number of complex roots 
organised into strings. These strings describe breathers in the sG language, and we 
shall not consider them. 

Changing the cycle, the DDV equation (17. 5p reads as follows 



(8.21 



log a(C) = log (^^^ - 2mud - ^{Cl^h) + *(^/^^^ 

- / /aC/01og(l + a(0)^ 



70 



where 



and 7o goes clockwise around M+. The first term in the right hand side is of order 
n. This implies that log a(C) is of order n. Standard analysis shows that 

(8.3) Re(loga(C)) > 0, Im < 0, 

Re(loga(C)) < 0, Im > . 



Hence for large n with exponential precision in n, only the lower part of the contour 
7o contributes and we come to the linear equation: 



(8.4) logaiO- I K{C/0\oga{0^ 



= log - 2mue - J2 *(c/a) + J2 "^(^z^^) 

We split log a(C) according to four terms in the right hand side: 



loga(C) = nF.a,(C) + FeiC) + F,{C) + F,(C) . 



The first term is model dependent which means that it depends on the inhomoge- 
neous parameters. But it does not contribute to physically relevant quantities. For 
example the S'-matrices for the XXZ model and for the sG model coincide. The con- 
tributions Fh{C) and -Fc(C) ^^e well known, so we shall not write the corresponding 
formulae. For F^IQ one finds immediately a ^-independent answer: 



" l—u 
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Let us consider the eigenvalue T{Q. Due to (18.31) with exponential precision in 
n, for < arg((^^) < vr we have 

logT(C) = log(a(C)) + mi,9 [ vl/(C/Orflog(l + a(0) 

Zm J 

7 

= log(a(C)) + muO - J2 ^(C/a) + E ^(^/^^) 

h c 

"2^/ ie"^^^'^^ log(a(0)rfe^ + O(e— ), 

where 

fCq-'^ - 1 
*(C)=log'^ ^ 



It is easy to see that 



c 



1 r d 



2 



SO we come to an important conclusion: in the limit n — t- oo with exponential 
precision in n, \ogT{C,) is independent of the twist 9. Finally a standard computation 
yields 



^•5) T(c)=nac(c)n^ 



where Tvac(C) is independent of 9. The complex zeros do not contribute. Returning 
to the picture fig. 4 we see that the dependence on a in the intermediate states is 
dropped in the limit n — )■ oo. The a-dependence remains in the quasi- local operators 

On the other hand it is clear that in this limit fig.4 represents just the form factor 
decomposition on the plane. Further, in the scaling limit the rotational symmetry 
occurs and we obtain the usual form factor decomposition for the sG model. 

9. Evaluation of w(C,0 in the presence of solitons 

We want to evaluate the function uj{(,C,) which corresponds to the right part of 
fig.4- In that case the Tl(C) is given by (18. 5p and T^{C) = ^vaclC)- Hence 

(9.1) ^>(c)=n^H^. 

where the number of holes is 2n. We shall use other variables 

c' = z, = 5, = . 

So, in our usual notation 

Later we shall use also r/" = H, ^-^ = X. 
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Usually the function u{(,C,) is found in two steps. First, we solve the equation 
(17. 7p for -Rdress and then substitute the result into fl7.9p to find uj{(,^). However, this 
procedure relies heavily on the assumption of simplicity of the spectrum of Tr(C)- 
In the case n — )■ oo this assumption is not true, the eigenvalue (18. 5 p corresponds to 
i^n) vectors in the sector of total spin 0. That is why the integral operator in (17. 7p 
is degenerate and -Rdress is not well-defined. To avoid this problem we proceed in 
the same way as [26] introducing the function G(C,0 which satisfies the equation 

(9.2) 

The shift e"*"*" in (19. 2 p refers to the poles of l/p{rf). The term 5'^ipQ{C,/^,a) as 
function of has poles at = ^'^ and = i'^(f- Let us choose such ^ that these 
poles do not lie on M+. Then the singular integral equation (19. 2p is well-defined 
being supplemented with the requirement that G(CjO regular for Q G M+. 
Once (^(C^O is found the function w(C,0 is defined by 

(9.3) a;(C,e) = <575.-A,-Vo(C/e,«) + T^ / 57^o(C/^, «)G(r/, 0" 



Writing down these equations we used (18.30 . 

The form factors were described in terms of pairings (£^"\ L*^"))q,. Now we want 
to solve similarly the equation (19. 2p . We formulate the result as Proposition, giving 
a necessary explanation for the precise definition of ^(CiO i^ ^^e proof. 

Proposition 9.1. Recall that 

n 

Mt\S^, ■■■,Sn) = n Sp U^Si - S',). 

p=l p>r 

Consider the rational function 

(9.4) i?a(5i, ■ ■ ■ , 5„) = -l zf^^,p^% n • • • ' Sn). 

^ ' p=\ p 

Then for any i^""^ the function 

(9.5) G{c,i-i^-^) = ) ; " 

solves the equation (19. 2p . 

Proof. First of all we have to explain how the pairing with the rational function 
Rz,x{SY"'^ is understood. For polynomials £(s), L{S) we defined the pairing as 

(£,L)„= J xicx)e^''^I^LiS)da + J xi^)e^''ni5)LiS)da, 

R-iO r 
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where 

(9.6) p(sq"2)£(s) = m(s) + a"2p(s)n(sq"^) - p(sq"^)n(s) , 

and m{s), n{5) are chosen from the requirement of convergence. Now for every Sp 
the rational function -R^^ poles at = Z'^. Admitting the existence of these 

poles we define the pairing (£("\i?^i')a by drawing the contours of integration as 
follows. 



Im A 




II] 



-2% 



Fig. 6: Contour of integration for the pairing. 

The presence of the poles at = Z"^ makes the pairing dependent on the choice of 
the functions m{s), n{s). 

We immediately see that there is a trouble with this definition. The point is 
that the polynomials m(s) and n{s) are defined from (19. 6p not uniquely. In spite of 
this arbitrariness the pairing was defined uniquely for polynomials L{S) because we 
could consider the difference between two solutions to (19.61) and cancel the integrals 
moving the contours. For rational functions with poles at S"^ = Z^ this is not true 
because moving the contours we pick up contributions from the poles {fig. 6). Thus 
the pairing depends on the choice of m{s) and n(s). Still we shall show that the 
equation (19. 2p is satisfied for any choice of m{s) and n(s) and, hence, the difference 
solves the homogeneous equation. 

Before doing the computation of the integral in (19. 2p . we give some remarks. First, 
it is very helpful for understanding of what is going on to check the last statement 
in the previous paragraph directly for n = 1 when we have only one integral over a. 
Second, we shall see that the arbitrariness in question does not concern the function 
a;(C, ^; •^'■"'■') which is the goal of our computation. The function G(^,^;£*^")) is an 
auxiliary object. It may depend on the choice of regularisation of the integrals, but 
the main object uj{C,^^] ^^"^) does not. This will be proved at the end of the section. 
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The last remark is regularity of G{(^,$,; i^"'^) at = ^j. It clearly follows from fig. 
6 that the integral has singularities at C = ij- They are simple poles, and P{Z) 
cancels them, so, altogether is regular at these points. 

We compute 



(9.7) ^ j {MQC/V,cy) - MQ~'C/V.cy)) R^H^xiSu - ■ ■ ,Sn 



drf 



1 
2wi 



\ 



i^o{C/V,(^)R-HxiSir-- ,Sn) 



( 



^o(C/^,«)<k5: 



Sn 



P{H) drf 
P{-H) r/2 



J 



X 



1) ■ ■ ■ ) 



- si 



Q2 _ C2 



Pi-X) 

MCS--", a)P{Sp) - AMCiSpQ)-", a)P{-S, 



PJ 5 



MCS-''e-^''"',a)P{Sp) - AMaSpQ)-''e-^^^'^,a)P{-S, 
where we used a trivial identity: 



PJ 5 



CTp G 



e r 



Im A 



4(logC+7a) 




-X-r^ 



4log^+7ti 



r;+2ra 



K- -X -X--^- ^'-X-^ 



-7a 



7tl 



4(iogc -™) 
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Fig. 1: Poles of the integrand in (19.71) . 

The position of the poles of the integrand is depicted in the complex (l/i^)log77- 
plane. 

The figure fig. 7 illustrates the position of poles. The poles in sequences headed 
by Pj and /3j + iri effectively occur when the integration in ap G Tj is performed. 

We are happy with the first two terms in the right hand side of (19. 7p . As for the 
last, containing the sum over p terms, we want to show that it vanishes when the 
integrals over ap are taken. Observing that in the summand for p the dependence 
on Sp before the curly bracket cancels out by one clearly sees the situation is 

similar to that considered in Proposition 13.41 for the Q-exact forms: both integrals 
over M and F reduce to residue at one pole like on fig. 2. The difference with the 
Q-exact forms is that ipo^C^Sp^ ,a) is not a 27ri-periodic function of ap. But on the 
other hand there is a difference between the integrands over M and F. It is easy to 
see that this difference is exactly such that the residues at the fat points on fig. 2 
cancel. This finishes the proof. □ 

The difference of two solutions is proportional to 



This gives a huge but finite number of linearly independent solutions to the homoge- 
neous equation. Indeed, each 4'^'' [i = 1, 2) is a linear combination of ii^i A ■ ■ ■ A ii^n- 
Every iij{s) should be split into mij{s) and njj(s). The degree of mij{s) is re- 
stricted by the requirement of convergence, and ?T-jj (s) is defined modulo p(sq^) (see 
(I3.22P ). Therefore, its degree is essentially bounded by 2n — 1. In addition Diff ((^, C.) 
is a polynomial of degree ra — 1 in X"^. 

Now we proceed to the computation of uj{C, 0- Using the function G{r], ^] we 
compute the integral in (19.31) . Since the dependence on rj and ^ is solely contained 
in -R^X' concentrate on the computation of the integral 



keeping in mind that the integration with respect to the variables Si, . . . , Sn is 
implied. 

To this end we want first of all to transform somewhat R^^\/ p{ri). Do the partial 
fractions 



(9.8) 




(9.9) 



i?i,"/^(^i,---,5„)Mr/) 
= ^m(")(5i,...,5„){ 



H H P{X) 



H-X H + XP{-X) 



X 



n 



X2 - S^. 



H 
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We call the last term in the bracket including the prefactor 2^Mq"'''(S'i, . . . , Sn) the 
sum term. Let us transform it by using the Q-exact forms: 

This is possible unless the point ^ log r] — ni coincides with one of the poles inside 
r. This does not happen if ^logrj is slightly above all I3j which can be harmlessly 
implied in the integral fl9.3p . So, our first goal is to compute the integral 



(9.10) 

1 

27n 



1 

2TTi 



V 



H 



A- 



H 



H-Sp H + SpQJ r/2 
C 



H dyf 



J 



A 



P{Z) 



Z - SpQ P{-Z) Z-Sp 



This computation is illustrated on the figure fig. 8 where the position of the poles 
and the effective poles is shown. The pole in the circle does not count because of 
the multiplier P{—Sp). We used A = e'^*" twice. 



Im 



-27t 



ilogC 




i(l0gC-7ti(l-V)) 



4(iogC-7ti) 



ap-27ii 



Fig. 8: Contour of integration for ( IQ.lOp . 



We use once again the Q-exact forms in order to transform 



A- 



Z — SpQ 



Z-S, 



-P{Sp 
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Thus, the sum term without the prefactor gives rise to 



P{-Z)P{-X) ^ S, 11 52 - 



z — s„ 



{P{S,)P{-Z) - P{-S,)P{Z)) 



Notice that the last part is a polynomial of Z and Sp. We divide it into the even 
and odd parts (see (15.41) for the definition of C{Z, S)): 

^^[p[Sp)P{-Z)-P{-Sp)P{Z))= -2uC{Z,Sp) 



+ 



Z ^ P{Sp)P{-Z)-P{-Sp)P{Z) _ P{-Sp)P{-Z)-P{Sp)P{ZU 



2 V 



Z-S^, 



j p Z -\- Sp 

For the second term the summation can be performed by the interpolation formula. 
We combine the result of this summation with the contribution to a;(C,0 coming 
from the first term in the right hand side of (19.91) : 



1 

2^ 



H 



H P{X) 1 dr] 



IH-X H + XP{-X) 



AuP{-Z)P{-X) 



/P{X)P{-Z) - P{-X)P{Z) P{-X)P{-Z) - P{X)P{Z)\ 

V z-x z+x / 



= -5^-5^-A;Vo(C/e,a)- 

We use the analytical continuation formula (I7.12p . Altogether we come to following 
nice expression 



(9.ii: 



(C,e;^^"^) 



Z,XJa 



where the polynomial i^^x given by 



(9.12) L^^USu---,Sn 



P{-Z)P{-X) 





X 

x^ 


C{Z,Si) ■ 
Si 

sf ■ 


C{Z, Sn 
Sn 
Sn 






C2n-1 



This final result does not depend on the arbitrariness of m(s), n(s) because L^^^^{Si, ■ ■ 
is a polynomial. 

10. Equivalence of BBS and BJMS fermions. 

Let us examine the result of the previous section. For finite n, the function C,) 
is associated with each eigenstate of the left Matsubara transfer matrix. In the limit 
n — i- oo, the eigenvalues are parametrized by a set of real numbers ■ ■ ■ ,/32n} 
and are (^")-fold degenerate. Correspondingly, for each fixed /3j's, one has to have 
the same number of the functions lu f*-"^). 
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In the process of solving the integral equation, however, we did not use any 
other condition on than that it is a skew-symmetric polynomial of degree at 
most 2?i — 1 in each variable. It enters the solution as a ratio, so the solutions of 
the integral equation are parametrized by points in the projective space P(/\" V) 
where V = ©^"q^Cs-'. Let us call ■ ■ ■ , /32n; ^'■"''1 the actual eigenvectors of the 
Matsubara transfer matrix for n — )■ oo. Then the corresponding ^'^"'^'s should be 
some special set of (^") points of the projective space. Let us see if one can further 
narrow down the possibilities. 

The main determinant formula fl7.4p tells that 

(10.1) det(a;(C.,0;^("))),^^.^,..., 





ib*(Ci)---b*(a: 


)c*(6) 


■■■c*(ei)g2"5(o)| 


vac) 




g2a5(0) 


vac) 



On the other hand, if is a pure wedge product £o A £i A ■ ■ ■ A of linear factors 
ii, then from the determinant formula fl9.lip - fl9.12p one deduces by a simple linear 
algebra that 

(10.2) 

det(a;(C e-^^"))) = (^^"^ ^ V^g(Zi) ■ ■ ■ ^g(Z,)xS(X,) ■ ■ ■ xS(^i)M^")). 

1 1 

Here, as usual Zj = , Xj = , and we used the definitions fl5.9p . flS.lOp . This 
hints at the following postulate. 

Postulate. The polynomials corresponding to the eigenvectors of the Matsub- 
ara transfer matrix for n — >■ oo belong to the Grassmannian Gr(n, V) C P(/\" V). 

We present below an argument in favour of this postulate, by invoking HLBA. 

For fixed n, the Bethe roots are either real or complex. According to [HT], in the 
limit n — )■ oo the real roots fill densely M+ with holes which correspond to solitons 
with rapidities /3j, and the positions of the complex roots are defined by the Bethe 
roots of the transfer-matrix constructed from the physical S-matrix. Namely, to 
the vector ■ ■ ■ , /32„; i^""^] there corresponds a Bethe vector in (C^)®^". Let us be 
more precise. The two soliton S-matrix is given by 

(10.3) SijiPi - Pj) = 5o(A - /3,■)^^,,(b,/b,■) , 

SM = exp I -I / sin(2fct./3)sinh((2t.-l)7rfc) 
/ccosh(7ri^A;) sinh(7r(l — z/)A;) 







and 
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We define 
(10.4) 



(c(t) -D(t)) ~ 'S'a,2n(t/b2n) " " " S'a,i(t/bi) . 



Then we are interested in constructing the eigen-covectors of the HLBA transfer- 
matrix 

T^^^^{i) = xA{i) + x^^D{t) , 

where x is such that \x\ = 1, it is defined by 6 from Section [HI The Bethe vectors 
in the weight zero sector are given by the algebraic Bethe Ansatz [19]. So, the 
correspondence between the eigen-covector in the n — oo hmit and the Bethe vector 
in HLBA reads as 

n 

(/3i,---,/32„;^(")|^(tl n^K-)' 

i=i 

where ( t I is the covector with all spins up, and u = (ui, . . . , , u„) satisfy the Bethe 
Ansatz equations. In this correspondence the eigen-covectors are parametrised by u. 
Suppose that /o„(/3i, ■ ■ ■ , /32n) is the form factor of Oq. Then, we have the equality 

n 

(10.5) ■ ■ ■ , P2n, ^('^^ |0«|vac) = ( t I n ^K-) ■ /(/^l' ■ ■ ■ ' '^2n) • 

i=i 

Let us compute the right hand side by using 02. 4p . The definition oiw'^^'"'''^^"{Pi, ■ ■ ■ 
is basically the same as in [7] , but it has to be transposed since we consider the ma- 
trix elements between the excited state and the vacuum and not vice versa: 



The overall multiplier and the integral transformation involved in the formulae for 
the form factors are independent of the partitions. So, it is easy to see that we have 
to compute 

fin) , 

1 " 

= E -— n v^(tinc'K) n c{bM)- 

{l,-,2n} .^t^ ' iel- j = l j- ^j= + 

=i-ui+ , 

The scalar product in the right hand side is the domain wall partition function given 
by the Izergin determinant [32], so, we face a difficult but clearly stated combina- 
torial problem. Surprisingly, for any set Ui, ■ ■ ■ ,u„ (not necessarily satisfying the 
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Bethe equations) we find that ^|"|(si, ■ ■ ■ Sn) belongs to the Grassmanian, namely, 

(10.6) 4:1(51, ■■■5.) 

n 

n (Ui -Ujq2) 

= C(q) ^'^"^ i{u},0 A i{u},l A ■ ■ ■ A £{«},„_l(Si, ■ ■ - Sn) , 

n n(u.-b,q^) 

i=li=l 

where c(q) is an irrelevant constant depending only on q, and 

n 2n 

n n (u. - b,q2) 

(10.7) £m,(5) = c(q)"^^^^^^^ 

n i^i - 

1 " 

X E n fb -b) nv^( tin n c'(b,)ii). 

=/-u/+ , 

Notice that 

^i-ui+,j = Cj,os""^^ H Cjjs" + rfj-is''"^ H h rfj> , 

where dj^k depend on the partition /"•" U and Cj^k do not. From this fact it is easy 
to see that the identity (110. 7p is a necessary condition for (110.61) to hold being a part 
of it. The rest of equations in the identity (110. 6p can be viewed as Pliicker relations. 

Now we see that our assumption about f'-"^ being in the Grassmannian fits com- 
pletely with the fact that the function u{C,^;i^'^^) must describe the ratio of the 
component of the form factor for the descendant to the one for the primary field in 
the basis of Bethe vectors. The formulae (110. ip and (110. 2p allow us to compute the 
form factors of any operator for homogeneous or inhomogeneous XXZ chain. Here 
we shall not consider the form factors for the lattice model and proceed directly to 
the sG case. 

As we have seen, for the description of the sG form factors, it is natural to slightly 
modify the operators ipo{C), XoiO by introducing a nontrivial block A in (I5.2ip . On 
the other hand, when we take the scaling limit from the inhomogeneous spin chain, 
the lattice fermions b*(C), c*{() are also modified by a Bogolubov transformation 
involving the function 

wo(c,o = -5c'^c"AcVo(c/e,«). 

If we are working with the asymptotic operators f3*{C), 7*(05 ^*{0^ 7*(0 this 
function contributes only to the pairings of /3*(C) with 7*(0, and of /3 (C) with 
7*(^). In these two cases we consider C — !■ oo, ^ — t- and C — ^ 0, ^ — oo respectively, 
and keep in the asymptotics (I7.18P of tpoiC/^i'^) o^ily the part going in (^/C)^ or 
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So, effectively a;o(C)0 is replaced either by 



oo 



or by 

oo 

Now it is easy to see that 

Due to this identity the effect of modifying h*{() and c*{() is exactly the same as 
the effect of modifying ip^iQ , X*oiO- 

Putting together all pieces of the puzzle we come to the main conclusion of this 
work: 
(10.8) 
^(*) 

/3*(MCi)---/3*(MCp)^*(Cp+i/M)---^*(a/M)r(?fc/M)---7*«,+i/M)7nM€.)---7*(M?^ 

So, the action of the BJMS fermions on local operators coincides with the action of 
the BBS fermions on towers. This exact identification of two things introduced orig- 
inally for completely different reasons is one more evidence of deep self-consistency 
of integrable two-dimensional quantum field theory. 

11. BBS CONSTRUCTION OF NULL VECTORS 

The main achievement of the paper [T7] consists in the fermionic description of 
null vectors. In this section we shall discuss this issue. Throughout this section we 
fix the number of the parameters Pj to be 2n. 

Up to now we were interested in generic a. However, as it is clear from the 
discussion of the regularised integrals, something special happens at the points of 
resonance, i.e., at the points where the pairing (13.171) has poles. In this paper we 
shall consider only the resonances occurring at a = ^-^m with m G Z>o. In the 
CFT language this corresponds to considering the degenerate fields which in the 
standard notation are 

$l,m+l = <^'i^,n- 

We shall use both of these symbols. The polynomials a"'i^pj^j+ depend on a only 
through a^, i.e., ^^-periodically. So, the form factors of the descendants O^.^i^ 
are expressible via the pairings ( , )a'- 



(11.1) (4"4,(5),lW(^))„^„^ = (-)-(4"J^,,(5),LW(^)n^: 

Having this in mind we shall consider all of them together. 



n 

j 
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In the present paper we shall consider only the case of the fields $i,2m which 
we call even. According to our logic all these primary fields and their Virasoro 
descendants should be considered as the fermionic descendants of 

$1^2 = ■ 

We denote this space by CKi^ = ©cez^^iizii c where CKi^ ^ denotes the sector of 
charge c. 

Similarly, the case of the fields $i,2m+i which we call odd is reduced to the case 
a = 0. While the point a = belongs to the fundamental domain (16. lip , a = 
is on its boundary, and this is a source of many complications. For that reason 
we decided to postpone the consideration of a = till a future work. Still some 
technical points explained in the next subsection will be common to both even and 
odd cases. 

11.1. Peculiar properties of form factors at = 1. Due to (111.11) . in order to 
treat the case = 1 it is sufficient to study 

(4"4„ivW(5)).L=„ 

for all Laurent polynomials N^'^\S). Let us do that assuming for the moment that 
the coefficients of N^^'^{S) are independent of a. 

Motivated by the formulae (I3.15p . (13.160 . (13.170 . let us introduce the following 
definition of residues for any Laurent polynomials N{S), ^(s): 

Res+[iV] = res5=oo(x+(5)5-2"iV(5)— ), tcs+li] = res,=oo (x+(s)5-"£(5)^) , 

Res_[iV] = Tess=o{x-{S)NiS)^), xtB.[i] = res,=o (x-(5)£(s)^) . 

Extending the definition of Res± as 
Res± [iViA---AiV,](5i,--- 

k 

= ^(-l)^'-iRes±[iV,-] (iVi A ■ ■ ■ A AT, A ■ ■ ■ A A^,) {Si, ■ ■ ■ , S^^i) , 
i=i 

we have 

(11.2) Res^ = 0, 

(11.3) Res±[i:>i[Z]] = 0, 

where Di[Z] = Da[Z]\^_^. Obviously the same relations hold for xts± defined 
similarly as above. Note that in the definition above, no multiple poles occur by 
antisymmetry. 

Then using (I3.15p . (I3.16P and (I3.17p . we can write the residue of the pairing as 
(11.4) 

Tes^=oi^^'\N^'^)^da = -i^((rc5+[£(^)], Res+[A^«])„ - (res_[£«], Res^iAT^^)]) J , 
for any i^'^^ and A^'-'^'^ whose coefficients are independent of a. 
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Consider the polynomials (.i-ui+,i = 0, 1, ■ ■ ■ , n — 1). We have 

^/-u/+,o(5) = (a"^ - a)q-2>^+(sq2) . 

The rest of them have the properties {\- < i <n — 1) 

(11.5) ^/-u7+,.(5) = aDj5'""] + 0(s""i) , 5 oo , 

^/-u7+,.(s) = (a"' - a)Ci{h) + 0(s) , s ^ , 

where Cj(b) are irrelevant constants depending only on bj's. These follow easily 
from f l4.3p . The formulae f lll.Sp clearly imply that 



te5±[£/-u/+,i(s)] 

Altogether we come to the formula 
11.6) 

lim(4"J,,, , iV(«))„ 



= 0, i = 1, ■ • ■ n — 1 

0=0 



= -2^^[{if-^l , Res+[iVW])„ + (-q)" H H ^'(^t-uA > Res_[iVW]). 
where 



is/- ie/+ 



a=0 



^/-U7+ — •^/-U/+,l A ■ ■ ■ A il~ul+,n-l ■ 

When a = 0, the monomial S"^"" cannot be reduced to lower degree because it drops 
from the exact form P{S) — AP{—S). So one has to relax the degree restriction 

dEED to 

(11.7) 0<deg5^iV(")(^i,--- ,^„) <2n 

An important point about a = is that for certain polynomials satisfying f lll.7p 
the pairing flll.6p vanishes. This happens if one of the following three conditions is 
met: 

(i) Vanishing of both residues. The pairing flll.6p vanishes if 

(11.8) Res+[A^(")] = Res_[iV(")] = . 

(ii) Exact form. The pairing vanishes if 

(11.9) iV(") = D[l] AiV("^i) , 
where 

D[1] = Da[1]U=, = P{S)-P{-S). 

(iii) Quantum Riemann bilinear identity. The pairing vanishes if 

(11.10) J^{n) ^ Qi2) ^ ^ 

where 

C^'\Si,S2)=C{Si,S2)-C{S2,Si). 
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The last property (iii) needs some comments. It is a consequence of the quantum 
Riemann bihnear identity [TB]. Consider the antisymmetric polynomial 



Suppose we find polynomials rj(s), Si(s) {i = 1, ■ ■ ■ ,n — 1) in the kernel of the 
operators 



res± : €[6]^"^^ ^ ©^^Y^C[b]s^' ^ C[b] 



such that 



n-l 



(11.11) C^^\5i,52) = iri{52)Sii3i) - r,(Si)Si(S2)) . 

Then a pairing in Ker res+ fl Ker res_ is defined by 

Ti O Tj = SiO Sj = 0, O Sj = 6ij . 

The sets {rj}"j/ and are called half-bases. Obviously there is an action of 

the symplectic group Sp{2n — 2). 

The quantum Riemann bilinear identity (see [16]) states 

(11.12) (mi A m2, C^^^)o = 27r?(mi o ms) . 

It is clear from (14. 3 P that for any partition {1, ■ ■ ■ , 2n} = U the following 
polynomials satisfy (111. lip for a = 1: 

r,(s) = (q^s)"-^ , Si{s) = ii-ui+A^) ■ 

So, 

(11.13) ii-ui+,i° ^i-ui+,j = , 2,j = 1, ■ ■ ■ ,n - 1 , 

which ensures that the pairing of ii-ui+,i A ii-ui+,j and C*^^^ vanishes. The pairing 
-UI+ A ^/-u/+ 1 and C*^^'' also vanishes due to 

Res+,5iC(2)(^i, ^2) = 0, Res_,s,C(2)(^,, ^2) = -^(P(52) - P(-S2)) 
and ffTTQl) . 

Before going further let us make one remark, which is supposed to be well known 
but probably worth being repeated. For generic a the polynomials •£/-□/+, A ■ ■ ■ A 
^/-u7+,n-i with different partitions of b/s span the space /\" C^" whose dimen- 
sion (^^) coincides with that of the weight subspace (C®^"')o. The latter can be 
interpreted in terms of the ?7(l)-symmetry of the sG model: soliton-antisoliton pro- 
vide a two-dimensional representation of U{1). When a = the polynomials 
•^7-u/+,i(s), ■ ■ ■ , ^/-u/+,n-i(s) which enter the form factor formulas are divisible by 
5, and res+ vanishes on them. So, they are linear combinations of 2n — 2 different 
monomials. Moreover, the form factors vanish if contains c*^^-* as a multiplier. 
Altogether, instead of /y'^C^", we have a smaller space /y"'"^ C^"~^ isomorphic to 
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the space of the maximal irreducible representation of Sp{2n — 2). The dimension 
of this space equals 

\n - I J ^ \ n - J ~ \ n J ^ \n - l) ' 

In the right hand side we have the Catalan number, the multiplicity of the singlet 
representation of the group SL{2) in C^^". This fact agrees with the quantum group 
reduction [12] which states that the form factors of ^i^^ are invariant under the 
action of the quantum group ?7q(s[2). Certainly, for generic q, this multiplicity is 
the same as in the classical case. 

11.2. Null vectors for the fields $i,2m- In the case a = and generic Q we 
can safely go from the odd representatives for L^"'^ to those of restricted degree and 
vice versa. So, following the remark at the end of Section |5] we consider the fermions 
defined by f l5.32p . which act on the space of towers with restricted degrees: 

0<degsp''\Su--- ,Sn)<2n-l. 

The identification with the previous subsection goes through 

n 

(11.14) N^''\S) = Y[Sj-L^''\S). 

As already noted, we have = (f^'^), iV("))o. 

First, observe that 

c{z, s) = ^z. ^(^) - P(-s) + o{Z') , 

2z/ S 

and 

C.{Z, X) = Z ■ ^ cot f (a - l)a2n{B)X-\P{X) + P{-X))l^^+OiZ') = OiZ') . 

It means that, under the identification flll.l4p . the operator ipl acts as wedge product 
by an exact form. Hence by (111.91) we find a set of null vectors: 

(11.15) V^*J{i^_i = 0. 

1/ ' 

We have 

n 

1 < degs^ (J[ S, ■ L(")(^i, ■ ■ ■ , Sn)) < 2n , 

which implies that 

Res_ ^n5,-L(")(5i,-- - ,S'„)j =0, 

n s, ■ L(")(^i, ■ ■ ■ , ^„) = n ■ (xi^^^^)^"'''"^ (5i, ■ ■ ■ , Sn-i) . 

i=i / i=i 
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This gives due to {xiY = another set of null vectors: 
(11.16) Xi^i^i^O. 

Now we turn to the Riemann bilinear identity. Introduce the operators 



We have for G _,_2 
(11.17) 

((Ceven + Ceven)-^) {Si, ■ ■ ■ , Sn-s) = ""C'cvcn ^ ^ ■ ■ ■ , Sn-s) ; 

where 

^even('5'l; 'S'2) = C{Si, 5*2)5';^ ^ — C{S2, 5'i)S'2 ^ . 

The only non-trivial part of this computation is to make sure that no cross-terms 
occur due to the Bogolubov transform (15.191) . (I5.20p . This follows from 

Since 

SiS2C^ll^{Si, S2) = —C^'^\Si, S2) , 
this together with (Ill.lOp and (I11.17P gives rise to new null-vectors 

(11.18) (Cevcn + Ceven) _2 - . 

In what follows we shall be interested in right-chiral null vectors. For the operator 
$1,2 it is easy to identify them: 

Xl ri^x}-Mt^ ^ 0, #(/+) = #(/-) + 1 
where we introduced the notation 

(11.19) =?/^*^ ■■■V;*^, x*i = X*a^ ■ ■ ■ X*a„ ^/=<r--^ap' X*i = xl^ ' ' ' xl,, 
for / = {ai, ■ ■ ■ , ttp}, fli < 02 < ■ • • < Op . 

We added the condition 1 ^ /~ in the second formula since we do not want to count 
the same null-vector twice. Notice that 

Consider the right-chiral descendants of $i,2m+2: 

rI.X*I-Xh,i„^)M^*^ , = #(/^) + m , 

where we recall that /odd('^) = {1, 3, ■ ■ ■ , 2m— 1}. At first glance it is not clear how 
to construct the right-chiral null vectors with Coven + Ccvon because X2j-i present in 
Ceven might spoil the product x*j ^^(m)- The solution was found in [T7]. Consider 
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According to the above considerations this can be transformed 

ernV;+X7-xL.M^^*^ =^ (-eeven)'"+V;+x/-X4,,MMW = 0. 

The latter identity is due to the fact that acting on x} ^^{m)^^*^ every Ceven can be 
replaced by the finite sum J2f=ii'2j+iX2j-i- 

Let us summarise. The following right chiral null vectors exist for $(-2^+1) i^- 

Xl rj^X^j-Xl^Mt^ , = + m + 1 , 

Introduce locally the notation 

iJfc = Span{^;+xI-MW | = A; , 1^/-}. 

It is easy to see that 

(11.20) iffc — )■ C^^pn-^-fc is an isomorphism . 

Using f lll.20p for = m + 1 we can combine the two kinds of the null- vectors into 

(11.21) eZtn ^;+xI-xL.M^o*^ ' #(^^) = -m-2, 

if 1 e /~ we get the first kind of null-vectors, and if 1 ^ J~ we get the second. 
For the left descendants the situation is even simpler because Ccvcn commutes with 

Ceven ^^X^V^L.M^O*^ #(/+) = #(/")- ^ - 2, 1 ^ J+ , 

which again can be put together as 

(11.22) Ceven X^- V^L.M ^0^^ ' = #(^") " " 2 

Actually, we have a symmetry between the left and the right chiral null vectors 
under the isomorphism 

(11-23) X*2j^l ^ ^2j-l > X2j-1 ^ i^2j-l , V'2j-1 ^ X*2j-1 , ^2i-l ^ X2j-1 • 

If 1 G /" in f lll.2ip and 1 G /"^ in flll.22p the identification simply follows from 
ffTT:20D with k = m + 1. If 1 ^ J" in inT2T\f and 1 ^ J+ in ffTr22|) the proof 
goes as follows. Denote by X the image of Ceven under f lll.23p . The operators 
(X, Ccvcn, [X, Ccvcn]) coustltute an s[2-triple, we have the isomorphism f lll.20p for 
k = m + 2, and CcvcnX is invertible on if_m-2- So, 

'YTJ 'YfP 'Y'N — l/om+2 rr pm+1 tt 

•^-n-m-2 — -^l^^even-^J '-'even -"ni+2 — '-even -"m-|-2 • 
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12. Comparison of null-vectors with CFT results 
If we specialise to a = ai,2m.+2 (m = 0, 1, 2, . . .) where 

ai,2m+2 = (2m + l)i^, 

the Verma module has a singular vector at level 2m + 2. Let W2m+2 be the 
singular vector, and W2m+2 the submodule generated by W2m+2- One expects that 
the null vectors in the sG model in the previous section should give rise to a fermionic 
description of the singular vector W2m+2 and the space y^2m+2- By using the results 
of [H , in this section we will check that it is indeed so (up to level 8 and modulo 
local integrals of motion). We plan to give some more details including the case of 
both chiralities in a separate publication. 
From fl6.6p one derives 

(12.1) /3>_2„.7}-+2™7L.M*a+2„.i^(0) = C^™(«)/3j+7}-7L.M*^nO), 

where #(J+) = #(J-) + m, C'M = (-)'"(™+^)/' 117=1 ^2,-i(«) C^(a), and for 
negative indices 

(12.2) f3\ = U2-ah,. 

Using (112. ip and the results of the previous section we obtain for the space 

oo 
fc=l 

the following fermionic basis 

{em,e.enr^'P*I+Yl-^^2rn+iy-^ , #(^^) = ^iH " 2m - 2 , 

where 



oo 



-m,even — ^ ^ /32j+l-2m72j-l+2m ) 



with the convention (112. 2p . 

For making comparison with CFT results, it is convenient to work with fermions 
I3a^^*, Ja^^* ^ normalised as 



= ^a(«)/3.^^^*, 7: = ^a(2 - a)7r^^*, 



■ ^ ' ' ' ' (^J!r(f + ^)' 

As it has been mentioned several times, at this moment we can identify the action 
of f3lj_i, 72j_i only in the quotient space 



k=l 
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As a basis in this space we shall take -P(l-2, 1-4; ■ ■ ■)^a where the P are monomials 
in the even Virasoro generators ordered lexicographically. Another basis in the same 
space is provided by monomials in the fermions, 

^ ai,--- ,ap,bi,--- ,bp Hai h'ap Ibi I bp ■ 

Identification of these monomials with the Virasoro basis (in V^") has been given 
up to level 6 in [1], and at level 8 in [28]. The general structure is as follows 

= (pr^(i_2, 1-4, ■ ■ ■ ) + daP:''\u, 1-4, ■■■))^a, 



where da = | v^(25 - c)(24A« + 1 - c), and the coefficients of P^"'', Pf"^ depend 
only on the central charge c and on the scaling dimension A^. The dependence on 
c is polynomial while the dependence on Aq, is generally rational, simple poles at 
certain negative integers may appear. 
The operator Cm,evcn reads now as 

P _ P'2i-i((2m - 1)^) ^j,^^ ^^p^ 

^2j-l(2 - (2m - 1)^) -,CFT CFT 

D2,_i(2 - (2m - iy-^)D,^_,^_,{2 - (2m - i) i_£)T'2.-i74^-2,-i , 

Actually, after the common multiplier r(z/)^^^v^ is extracted, the remaining combi- 
nations of F-functions in the right hand side collapse to rational functions of u. 
The following are the fermionic null vectors in W2m+2 up to level 8: 



^5,3, ^1,3,3,1 



level 2 










level 4 


*3,1 








level 6 


^5,1, 


^1,5 - 


4(1/2 
1/2- 


-4) 
■16 


level 8 


^7,1, 


*1,7- 


9(1/2 
1/2- 


-4) 
■36 


level 4 


^1,3 








level 6 


^3,3 








level 8 


*5,3, 


^1,3,3,1 






level 6 


^1,5 








level 8 


^3,5 








level 8 


^1,7 









In every module the null-vectors of lowest possible degree are singular vectors. 
Modulo the integrals of motion, the simplest singular vector is 

W2 = (l-2 - Y^^-l)'^h2 = 1-2$1,2 , 

which agrees with the result [1] 

^^1,1 = 1_2$1,2. 

The next singular vector 104^ is less trivial. 

"^"(2^-^- 3(1-.) '-^n'^- 
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On the other hand we have for generic a [T] 

(12.3) vl>, 3 ^ [-A\ + ^ ^U) 

Substituting a = ai^ we indeed find that \l/i^3 = W4. 

The formulae for wq and on one hand, and for \l/i^5, 7 [H |28] on the other, 
are much more comphcated. Still, substituting a = ai^, a = ai^g into the latter we 
find perfect agreement. 

Consider the first null- vector which is not the singular vector in W2: 

It is easy to check that (1/2)1_2W^2 coincides with the formula from [1] 

c-lQ + 3da, \^ 
*""(2'- + -9 H*"' 

specialised to a = ai^2- 

Proceeding in the same way we have checked that all the vectors from the above 
table indeed coincide with the Virasoro null- vectors. In every particular case we have 
to solve an overdetermined system of linear equations, and eventually find solutions. 
We consider this fact as a strong support to the statement that the BJMS fermions 
when they are identified with the BBS fermions create the null vectors in the Verma 
modules. 
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